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In  a previous  report  we  studied  rational  approximations  to  the 
incomplete  gamma  function.  These  were  based  on  the  asymptotic  expansion  of 
the  latter  function.  In  the  present  study,  a similar  analysis  is  done  for  the 
same  function  based  on  one  of  its  ascending  series  representations.  The  work 
is  related  to  ^tfe^  developments  -in  £2]  and  {3^ , but  there  are  numerous  improve- 
ments and  seme  new  features.  


I.  THE  INCCMPIETE  GAMMA  FUNCTION 


We  consider 


i 

? 


In  (1.3)  the  path  of  integration  lies  in  the  z plane  cut  along  the  negative 
real  axis  and  is  the  ray  T)  exp(iu)  , T)— except  for  an  initial  finite 
path.  If  z ^ 0 , |o|  < n/2  , the  integral  in  (1.3)  exists  without  the  re- 
striction on  i>  . 

The  connection  of  these  functions  to  the  confluent  hypergecmetric 
functions  is  given  by 


§(a,c;z)  = jFj^ajcjz) 


(1.5) 


r(z,u)  = zue"zt(l,l+v;z)  = e"zf (l-v,l-vjz)  , 


(1-6) 


ijr(a,c;z)  = — dil£i  $(a,c;z)  + r^~-^  z1"c$(a-c+l,2-cjz ) 
T(a-c+l)  f(a) 


(1.7) 


The  statement  (1.4)  is 


y(.,.)  - i 7-f-  ■ i • 

k=0  (u^k+l  k=0  k.(w+k) 


(1.8) 


In  the  above  formulas,  standard  generalized  hypergeometric  notation  is  used. 
Thus 


F ("al'a2'“*'ap|zN)  = 
^ '>^1*^2,  • • .,bQ  * / 


“ pFq(al> *2>  • • * • •>tq>z  ) 


a>  IT  (ai  )fczk 

Z ^ , (a)k  = 

*r0  * k r(a) 

rr  (WJ 

i=l 


(1.9) 


For  further  discussion  of  the  hypergeometric  functions,  confluent  hypergeo- 
metric functions  and  the  incomplete  gamma  function,  see  [4,  Chs.  2,4, 6j, 

[5,  Ch.  9],  [6],  [7],  and  [a] . 


f; . j . •. ; t > f‘ 


Also, 


P (z,u)  = (~)n+1.e'Z  f (z-t  )ntn+vetdt 

Zu(u+l)n  “o 


, »n+l  . «z  2n+l 

= (-)  n'e  z 


(u+l); 


2n+l 


1F1(n+u+lj2n+u+2;z ) , 


Rn(z,«) 


(- )n+1n!r(u+l)r(n+u+l)z2n+11F1(n+lj2n+u+2;-z ) 
r(2n+v+l)r(2n+v+2 ) ^F^(-n;-2n-w jz ) 


(2.4) 


(2.5) 
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Both  A (z,  d)  and  Bn(z, u ) satisfy  the  same  recurrence  equation 


(2n+u+l)(2n+u+2) 


(2n+u )(2n+u+2 ) 


+ 7 T7 rAnl(z,v)  • (2.6) 

(2n+u )(2n+u+l)  n"^ 


Me  list  some  other  useful  relations : 


(2n+v)DBn(z,u)  - nBn(z,v)  + nBn_x(z,u)  = 0 , (2.7) 

(2n+u  )DAn(z, v )+(n+v  )An(z, u )+nzAn_L(z, v V z_1u(2n+u  )[An(z, u )-Bn(z, v )] = 0 ,(2.8) 


[zD2-(z+2n+u)D+n]Bn(z,v)  = 0 , D = — 


(2.9) 


Consider 


Me  now  develop  a useful  estimate  of  the  remainder  Rn(z,u)  . 


y = e“z/21F1(a, c;z)  , zy"  + cy'  + (K-z/4)y  = 0 , K = c/2  - a . (2.10) 


Assume 


, a (z)  = 1 , u = i(c-l) 


k=0  u* 


Then  the  a^s  can  be  generated  by  the  expression 


ak+i(z)  = -iz  ak(z ) + i f (t/4-K)ak(t)dt 


(2.11) 


(2.12) 


Clearly  the  series  (2.11)  is  absolutely  convergent  for  all  z as  it  can  be 
found  by  multiplication  of  the  series  ^F1(ajc;z)  by  the  series  for  e"z/2  . 
However,  if  z is  fixed,  c is  large  and  K « c , then  (2.1l)  is  a useful 
representation  of  y for  large  c . The  expressions  (2. 6)- (2. 8)  can  also  be 
deduced  from  some  general  results  of  Olver,  see  the  discussion  in  [8,  p.  76]. 

Apply  (2.10)-(2.13 ) to  the  confluent  functions  in  (2.5).  Then,  with 
the  aid  of  the  duplication  formula  for  gamma  functions,  we  get 


(- )n+*n!nr(v+l )r(n+v+l)z^n+^e”z  T ^ 


»vi)[r  (n.  s|i)  r („♦  I +1)]  z 

k-0 


nn!r(u+l)r(n+,.+l)z2n  ^e  z 


24n+2,,(2n+u+1)jr  (n+  T (n+  | +l)j2 

X|l+0(l/n3)}  , (2 

wnere  u = n + ^(v+l)  > K = -u/2  and  the  a^'s  are  given  by  (2.12).  It 
follows  that  for  z and  v fixed, 

lim  Rn(z,v)  = 0 (2 


Since  r(n+a)/r(n+b)  = na“^  [l+0(l/n)J 


/ xn+1  . . * 2n+l  -z 

R (z,  v ) - Cli nHulXlz e_ 

24n+2u+lnu(n, f 


, it  is  convenient  to  write 
{l+0(l/n)j  . (2.17) 


. . Concerning  the  use  of  (2.l),  it  is  important  to  know  the  nature  of 

z^y^  , i = 1,2, . . .,n  , the  zeros  of  -Bn(z,u)  . It  is  clear  from  (2.9)  that 
the  zeros  are  simple.  If  u = 0 , Bn(z, u)  is  essentially  the  modified  Bessel 
function  of  half-odd  order  (see  (4.2)).  In  this  instance  the  zeros  can  be 
deduced  from  the  work  of  01ver[9].  See  also  Grosswald  [l2] . The  zeros  lie  in 
the  left  half  plane  and  are  always  complex .except  when  n is  odd  in  which 

case  there  is  a single  real  root.  If  z^u ' is  the  magnitude  of  the  smallest 

root(s)  of  Bn(z,v)  , then  from  the  work  of  Olver  [9]  z^°) ~1.32548n  . Here 
2(t2-i)2  = 1.32548  where  t0  is  the  real  zero  of  t = coth  t . Thomson  [lo] 
and  Kublanovskaia  and  Smirnova  [llj  have  tabulated  gz^  q , the  former  for 
n = l(l)9  to  4d  , and  the  latter  for  n = l(l)21  to  5d  . * Salzer  [l3]  has 
tabulated  z^“l)  for  n = l(l)l6  to  15d  . We  have  prepared  some  tables  of 

for  u = -3/2, l/2,l  . For  all  the  u values  mentioned  above,  the  rr  jts 

life  in  the  left  half  plane.  If  n is  odd,  there  is  a single  real  root,  other- 

wise the  roots  are  complex.  On  the  basis  of  the  data,  we  conjecture  that 
z^u ' ~<1.32548n  + u + 1 - l/n  . An  alternative  conjecture  is  that 

Z2n  ~ ?n  2 (3n+u+l)  , z2n+l  ~ V2  nn2  + u + 2 . Thus  as  n increases  , the 
magnitude  of  the  smallest  zero(s)  of  Bn(z,u)  increases  linearly  with  n . 

We  see  from  (2.17)  that  to  achieve  high  accuracy,  n must  be  considerably 


larger  than 


and  so  the  value  of  z 


(u) 


is  not  critical. 


Another  rational  approximation  for  y(zein, v)  is 


uz_ue'z’luTTy(ze?-n,u)  = Sn(z,u)  + TQ(z,u) 


Dn(z,u)  Dn(z,v) 


(2.18) 


i 


C (z,v)  - - * T zn  F /'-n+l+k, n+u+l+k, II  1\  (2.19 

z k =0  (v)k+10t+D!  3 k+2  \ z) 


and  Dn(z,  v)  is  zn  times  the  jFj_  in  (2.19)  vith  k+1  set  equal  to  zero 
Thus  Dn(z,u)  is  Bn(z,v)  if  in  the  latter  we  replace  v by  v-1  . Also 


.n  -v  -z 


Qn(z,v)  = LI  e_„  J (z.t  ^n+^-lgt^ 


/ %n  , -z  2n 

(u+1)on 


^F^  ( n+  u ;2n+  u+1  j z ) 


(2.20 


Both  Cn(z,v)  and  D^z^v)  satisfy  (2.6)  if  v is  replaced  by 
v-1  . Also,  Dn(z,v)  satisfies  (2.7)  and  (2.9)  with  v replaced  by  v-1  • 
The  relation  analogous  to  (2.8)  reads 


(2n+v-l)DCn(z, v)  + (n+v-l)Cn(z, v)  + nzCn_1(z, v) 


+ z_1v(2n+v-l)[cn(z,v)-Dn(z,  v)]  = 0 


After  the  manner  of  deriving  (2.15)  and  (2.17)  we  have 


(2.21 


f-  )nTm.*r(v+l  )r(n+v )z2ne”ze2a^z 


Tn(z,v)  = 


,4n+2v-2 


(2n+» )[r  (n+  |)  r (n+  ^ )]! 


[l+0(l/n3 ) } 


)nnl'{u+l  Jz^ne~ze^a^~<Z  >/u*  (1+0(j/„)}  , 

24n+2v-lnv-l(n:)2  *•  1 


(2.22 


where  u*  = n+v/2  , a£(z)  = z(z+4v-4)/l6  , whence  for  z and  u fixed, 
lim  Tn(z,v)  = 0 . 


h 


If  z and  u are  positive  and  fixed,  then  for  a given  n , it  is 
clear  that  Rn(z,y)  and  Tn(z,y)  are  of  opposite  sign.  This  leads  to  the 
inequalities 


Cn(z>v)  -v  -z  r vl  t.  > An(z,v ) , » 

D»(2,»)  { dt  < ’ ( 2 ’ 


where  > or  < sign  is  chosen  according  as  n is  odd  or  even,  respectively. 
For  example. 


l_u+l)(u+2)-z  < uz-ue>z  r tu-letdt  < 
(u+l)(u+2)+z 


w+1 


U+l+Z 


, z > 0 , u > 0 , (2.24 ) 


— < e_Z  < — , z > 0 , 

2+z  z+1  ' ' 


(2.25) 


15-4z2  -1  -z£ 

< z e 


15+6z£ 


/ 


dt  < 


2z2+3 


z > 0 


(2.26) 


In  the  last  three  equations,  equality  prevails  if  z— >0 


At  this  juncture,  we  present  a summary  of  the  material  in  the  re- 
maining sections  of  the  report.  If  v = \ , (2.1)  yields  approximations  for 
the  error  function  and  related  integrals  which  is  the  subject  of  Section  III. 
Approximations  for  the  exponential  function  and  the  circular  functions  which 
come  from  (2.1)  when  u = 0 are  discussed  in  Section  IV.  Tables  of  the  poly- 
nominals  in  the  approximations  of  the  functions  discussed  in  Sections  III  and 
IV,  and  related  data,  are  found  in  Appendix  A.  If  y * 0 , (1.3)  is  the  ex- 
ponential integral  and  in  [l],  we  developed  rational  approximations  for  this 
function  based  on  its  asymptotic  representation.  Though  r(z, 0)  is  defined  by 
(1.3),  clearly  y(z, 0)  is  not.  However,  we  can  give  an  ascending  series  rep- 
resentation for  a function  closely  related  to  T(z,0).  We  show  in  Section  V 
how  rational  approximations  to  this  related  function  may  be  derived,  and  tables 
of  the  polynomials  in  these  approximations  and  related  data  are  given  in 
Appendix  B.  FORTRAN  codes  for  computing  rational  approximations  to  the  in- 
complete gamma  function  and  its  special  cases  are  presented  in  Appendices  C, 

D and  E. 
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III.  ERROR  FUNCTION  AMD  RELATED  INTEGRALS 


We  list  below  formulae  useful  for  the  approximation  of  the  error 
functions  and  the  Fresnel  integrals.  These  are  based  on  (2.1)  for  i>  = a . 

It  is  clear  from  (2.1)  and  (2.18)  that  Vn(z,£)  and  Rn(z,^)  can  be  replaced 
by  Sn(z,^)  and  Tn(z,-|)  , respectively. 


z x 

Y (2>s)  - J t“2e-tdt  = 2z5e"z  f Vn(ze’in,5)  + Rn(ze-in,^)}  . (3.1) 

0 J 
z 2 j 

Erf(z)  =^y(z 2,-£)  = J dt  = (n/4 )2  - Erfc(z)  , (3.2) 


Erfc(z)  = J1  e-t  dt  , 
z 


(3.3) 


Z 

Erf(z)  = ye"t2dt  = ze’z2  {vn(z2e_lrr,i)+Rn(z2e":lTT>i)]  , (3.4) 


Erfi(z)  = f et2dt  = - iErf(iz)  , 
0 


13n  i3n  z 

Y(ze  2 >h)  = e4 

'0 


in  in 


J'  t"2e4tdt  = 2z2eiz"t’i3TV,4^n(ze^  ,^)+Rn(ze2  ,^)  ^ , 

z 

C(z ) + iS(z ) = (2n)"2  f t“^eitdt  . 


(3.5) 

(3.6) 

(3.7) 


The  exact  coefficients  of  polynomials  simply  related  to  AQ(z,^)  , 

> cn(z>£)  and  Dn(z>£)  are  given  in  Appendix  A £or  n = 0(l)l0. 

For  n * 2(1 )10  , we  also  record  values  of  R*(z,^)  * 1 2z2ezRn(z,^)  | , see 
(2.15),  for  z = re4®  , r = 1(1)10  ,9=0,  n/2  , n . Error  tables  for  other 
values  of  u and  for  Tn(z,u)  are  easily  derived  from  the  latter  table.  See 
the  discussion  in  Appendix  A.  See  Appendix  C for  FORTRAN  codes. 
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To  illustrate  the  utility  of  (2.15),  ve  present  the  numerics  below. 
The  tables  give  e_i'T/2Y(zein,^)  , its  rational  approximation  (2.1),  the  exact 
error,  and  the  approximate  error  according  to  (2.15).  The  data  are  developed 
for  z = 2e19  , n = 4 . 


f t"?etdt 
0 

6 Exact 


0 6.6876855 

n/2  3.3756999  x 10-1  + 2.3328174i 

n 1.0324115  x 10-8  + 1.6 91806 7i 


C Exact  Error 

0 -1.47  x 10-5 

. i/2  8.31  x 10"6  - 6.70i  x 10’6 

TT  6.16  x 10-12  + 1.201  x 10-5 


2zkzVn(z,i) 

6.6877002 

3.3756168  x 10"1  + 2.33282411 
1.0317960  x 10'8  + 1.6917947i 

2z2ezRn(z,i) 

Approximate 

-1.47  x 10”5 

8.35  x 10-6  - 6 . 74i  x 10-6 
1.05  x 10”11  + 1.19i  x 10-5 


IV.  THE  EXPONENTIAL  FUNCTION 


If  u — >0  , (2.1)  becomes 

e'2-s£r  + R”(z'0)  ' (4-1) 

Gn{z)  * znpF0{-n,7i*l;-i)  - (2/,T)5znKn+Jz/2)  , (4.2) 

eZ^i+^(z/£ ) n[(i)n]2 


-Li 


n+l  2n+l„-z 


ttz  ~e 
4n+1(nJ  f 


{l+0(l/n)} 


(4.3) 


Here  In(z)  and.  Kj^z)  are  the  familiar  notations  for  the  modified  Bessel 
functions . 

It  is  of  interest  to  show  how  (4.1)  can  he  used  to  compute  the 
exponential  and  circular  functions  in  an  efficient  manner.  Let 


Gn(z)  = Mn(z2)  + zNn(z2) 


where  M^z2)  and  Nn(z2)  are  even  polynomials  in  z . Clearly 


(4.4) 


rz  - Gn(~z ) _ M^Z2  )**zNq(z2  ) 
n Gn(z ) M^z2  )+zNn(z2 ) 


(4.5) 


and  one  readily  verifies  that  Gn(z)  , Mjj(z2)  and  Nn(z2)  all  satisfy  the 
recurrence  formula 


Gn+i(z)  = 2(2n+l)Gn(z)  + z2On-1(z)  . 


(4.6) 


Let  the  approximation  to  the  circular  functions  be  denoted  by 


e"iz  = cosn  z - i sinn  z 


(4.7) 


Then 


COSn  z = 


-2a^2,  .1^  . ■ ?.4f± 

Wn(z2)  Wn(z2) 


un(z2)  = [^(-z2)]  -*2[*t(“Z2)]  , Vn(z2)  = ^(-z2  )Nn(-z2)  , 


Wn(z2)  = [Mn(-z2)]2  + z2[Nn(_z2)]2 


(4.6) 
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It  can  be  shown  that*  UQ(z2)  , Vn(z2)  and  Wn(z2)  all  satisfy  the  recur- 
rence formula 

(2n-l)Un+1(z2)  = [-z2+4(4n2-l)j  [ (2n+l)Un(z2 )-z2(2n-l)UQ.1(z2 )] 

+ z6(2n+l)Un_2(z2)  . (4.9 

The  exact  coefficients  of  the  polynomials  G (z)  , Un(z2)  , Vn(z2)  and 
V/n(z2)  are  given  in  Appendix  A for  n = 1(1)10.  See  also  the  introduction 
to  Table  A. Ill  in  Appendix  A for  the  evaluation  of  Rn(z,0)  . Appendix  D 
gives  FORTRAN  codes  for  the  computation  of  e"z  , cosn  z , sir^  z and 
tann  z for  z real  only.  We  conclude  this  section  with  an  example  to  in- 
dicate the  effectiveness  of  (4.l)-(4.3). 

If  n * 4 and  z = 1 , (4.5)  gives  the  value  e^1  = 0.36787  94561 
The  error  is  -1.5  x 10"8  whereas  the  last  of  (4.3)  gives  the  value 
-1.6  x 10'8  . For  n = 4 and  z = i , (4.5)  yields  e*1  = 0.54030  23380  + 
0.8414709642i  . The  true  error  is  -3.2  x 10“8  - i2.1  x 10“8  as  compared 
with  the  estimated  error  from  the  last  of  (4.3),  -3.4  x 10"8  - i2.2  x 10“8. 

z 

V.  RATIONAL  APPROXIMATIONS  FOR  E(z)  = z"1  f t-1(l-e‘t)dt 


In  this  section  we  develop  some  rational  approximations  for  E(z) 
It  is  of  interest  to  first  list  some  functions  related  to  E(z)  . We  have 

z 

E(z)  = z"1  J t“1(l-e"t)dt  = z-1[-Ei(-z)+Y+ln  z]  , -Ei(-z)  = T(z,0)  , (5. 

where  y is  Euler's  constant. 


* We  are  indebted  to  Mr.  Jet  Wimp  for  proof  of  this  statement  and  other 
helpful  suggestions. 


(5.2) 


j[E(zein/2)+E(ze~iTT/2)]  = z‘1Si(z)  = z'1  J t-1sin  t dt 

z 

= z-1[tt/2  + si(z)J  , si(z ) = J t_1sin  t dt  , 

oo 


2i[E(zeiTT/2)-E(ze"iTT/2)]  = 4z_2H(z)  = 4z‘2  J'  t"1(l-cos  t)dt 


_ z 

= 4z’2[ci(z)-y-1h  zj  , Ci(z)  = / t-1cos  t dt  . (5.3) 


Approximations  for  E(z)  can  be  derived  on  the  basis  of  the  results 
given  in  Section  II.  Clearly  from  (l.l)  and  (2.1), 


z 

zE(z)  = lim  J tu-1(l-e_t)dt 
u — *0  0 


lim  - y(z»v)1 
v— ♦ Cr  J 


be  approximated  by 


(5.6) 


where  f^(z)  and  gg(z)  are  polynomials  in  z of  degree  p and  q , 
respectively.  If 

gq(z)E(z)  - fp(z)  = zP+^+1h(z)  , h(0)  t 0 , (5.7) 

then  (5.6)  is  the  Pade  approximant  of  E(z)  . Thus  (2.1)  is  of  the  type  (5.6) 
where  p = q = n , and  (2.18)  is  also  of  the  same  type  with  q = p+1  = n . 

The  approximation  (5.6)  is  called  the  main  diagonal  Pade'  approximation  if 
p = q = n . 

The  numerator  and  denominator  polynomials  in  the  Pade/  approximants 
to  transcendental  functions  are  known  in  closed  form, as  for  the  functions  in 
Section  II, only  in  very  few  cases.  However,  the  Pade'  approximant  to  a Taylor 
series  expansion  can  always  be  found  by  solving  systems  of  linear  equations. 
Thus,  if 

f(z)  = 21  fkzk  ' «(*)  = 21  8kzk  > (5.8) 

k=0  k=0 


then 

r 


p+q+l+k 


•j6p+q+l+k-j 


and  h.  must  be  determined  anew  for  each  choice 


As  remarked  previously,  it  is  inconvenient  to  use  (5.4)  as  a 
rational  approximation  for  E(z).  To  circumvent  this  difficulty,  ve  have 
computed  the  coefficients  of  the  main  diagonal  Pade'  approximants  of  the 
functions  E(z)  for  n = 0(l)l0  and  z“^Si(z)  and  4z"^H  (z)  for 
n = 0(2 )10  . These  are  tabulated  in  Appendix  B. 

We  also  include  tables  of  the  absolute  values  of  the  errors  incurred 
by  using  the  Pade'  approximation  of  E(z)  for  n = 2(1)10  , r = l(l)l0  , and 
9 = 0,tt/2,tt  , where  z = re1-®  , and  the  Pade'  approximants  of  Si(z)  and 
H(z)  for  z real,  z = l(l)lO  and  for  n = 4(2)10  . These  tables  may  be 
used  as  a guide  in  selecting  the  order  of  approximation  necessary  to  obtain  a 
desired  accuracy.  For  example,  if  six  decimal  accuracy  i6  desired  for  E(z ) 
for  z = Zzel'n/l  > i.e.,  I error  I < 0.5  x 10"®  , interpolation  of  Table  B.II 


indicates  that  a third  order  approximation  should  be  sufficient.  The  third 
order  approximation  gives  0.76507  22371  - 0.15899  838671  and  the  true  value 
is  0.76507  22539  - 0.15899  862561  . Thus  | error  I < 2.41  x 10'7  . 


E(iz)  = z_1Si(z)  - iz“^H(z) 


(5.10) 


Select  n . It  is  readily  deduced  from  the  error  tables  that  the  Pade'  values 
for  E(iz)  are  better  than  the  values  deduced  from  (5.10)  by  using  the  Pade' 
approximants  for  S(z)  and  H(z)  for  z real  . Thus,  if  both  Si(z)  and 
H(z)  are  needed,  it  is  better  to  use  the  Pade  for  E(iz)  . However,  if  only 
H(z)  , say,  is  needed,  it  is  more  economical  to  use  the  Pade  for  H(z)  . 


An  examination  of  the  zeros  of  the  denominators  of  the  Pade  approxi- 
mations of  E(z),  Si(z)  and  H(z)  indicates  that  the  magnitudes  of  the 
smallest  zeros  of  the  denominator  increase  linearly  with  n . Since  n must 
be  significantly  larger  than  the  argument  to  attain  good  accuracy,  the  loca- 
tion of  these  zeros  is  not  important. 

As  noted  above,  the  Pade'  approximants  for  E(z ) are  not  known  in 
closed  form,  and  it  is  necessary  to  solve  systems  of  linear  equations  for  each 
selection  of  p and  q (see  (5.8)).  We  now  give  another  rational  approxima- 
tion of  E(z)  which  is  of  the  form  (5.8)  with  p = q = n . In  over-all 
accuracy,  it  is  somewhat  inferior  to  the  corresponding  Pade  approximant,  see 
Tables  B.II  and  B.VIII.  However,  it  has  the  desirable  advantage  that  the 
numerator  and  denominator  polynomials  can  be  computed  by  recurrence  formulas. 
Following  [l4j , it  can  be  shewn  that 


1 2 (n+1 ) J e(n+l)2(2n-l)  J 


and  so  for  z fixed. 


lim  en(z)  = 0 , (5.15) 

n— 


whence  the  rational  approximants  converge.  It  will  also  be  demonstrated  else 
where  that  both  9q(z)  and  fQ(z ) satisfy  the  recurrence  formula 


fn(z)  = (B1z+A1)fn.1(z)  + (B2Z+A2  )fQ_g(z ) + Ajfn_j(z)  > 


A,  = - A,  = A—  = X , 

1 3 n(2n-3)  ' 2 


Bl  = , and  B2  = -2.(?p-1  )(ar?i 


(5.16) 


To  illustrate  the  power  of  (5.14),  take  z = \ and  n = 4 . The 
true  value  (from  (5.12))  is  1101  and  the  value  using  (5.14)  is  1098.5. 

From  the  preceding  development  and  the  last  remark,  it  is  obvious 
that  all  that  is  needed  for  a useful  expression  for  the  error  ien(z)  | is  an 
approximation  of  Fn(z)  (see  (5.1l))  for  large  n . This  is  not  available 
at  the  present  time.  If  z is  real  and  positive,  we  can  show  that 

(5.17) 


This  bound,  however,  is  very  conservative. 

The  coefficients  of  the  polynomials  ^(z)  and  fn(z)  613:6  6iven  in 
Appendix  B for  n = 0(l)l0  . Also  presented  are  values  of  |en(z)  I for 
n = 3(l)l0  and  for  z = re1-^  where  r = 1(1)10  , 9 = 0,tt/2,tt  . 

The  approximation  (5.11)  has  the  same  properties  as  the  Pade<  ap- 
proximation in  that  the  magnitude  of  the  smallest  zeros  of  the  denominator 
polynomials  increase  linearly  with  n , the  order  of  approximation.  Again 
since  the  order  of  approximation  must  be  significantly  larger  than  the  argu- 
ment, the  location  of  zeros  of  the  denominator  polynomials  is  not  critical. 


» 


I £ 

5 fat 

fn(z) 
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APPENDIX  A 


COEFFICIENTS  FCR  RATIONAL  APPROXIMATIONS  TO  THE  ERROR  FUNCTION 
EXPONENTIAL  FUNCTION  AND  CIRCULAR  FUNCTIONS 


We  first  show  how  the  exact  coefficients  in  the  approximations  (2.1) 
and  (2.18)  can  be  generated  when  v is  rational.  Suppose  v = p/q  , q ^ 0 
and  p,q  are  co-prime  integers. 


A*(z,u)  = q2n 


b£(z,v)  = q2n 


r(»+i) 


r(v+i) 


Vz>”)  = Z an,kzk  > 

k=0 


Bn(z,u)  = E b kzk  , 
k=0 


(A.l) 


where  An(z,v)  and  Bn(z,v)  are  defined  in  (2.2)  and  (2.3),  respectively. 
Note  that 


An(z, w ) _ An(z, u ) 
B*(z,v)  BQ(z,v) 


(A. 2) 


and  the  transformation  (A.l)  insures  that  the  coefficients  an  k and  k 
are  integers  if  p and  q are  co-prime  integers.  If  v = 0 , set  p = 0 
and  q = 1 . Frcm  (2.6)  it  is  seen  that 


A*(z,u)  or  B*(z,v) 


Aq(z,u)  = 1 , A-^z,!;)  = (p+q)(p+2q)  - qdz  , 
b*(z,u)  = i , b£(z,v)  = (p+q)(p+2q)  + q(p+q)z 

Using  (A. 3)  and  (A.l),  we  get  the  recurrence  formula 


l,k  = (2nq+p+q  )(2nq+p+2q  )ln^  ^ Xn,k-1 


■>n(2nq+p+2q ) (nq+i 


(2nq+p ) 


n-1, k-2  * 


xn,k  = an,k  or  bn,k  i k = 0,1,2,... ,n 
and  ]j  = 0 if  k<0  or  k>n 

The  initial  values  are 

ao,o  = 1 ' 

a^0  = (p+q)(p+2q)  , ai;i  = -<i2  , 

bo,o  ^ * 

t»1)0  = (p+q)(p+2q)  , bi,i  = *i(p+(i) 


Using  a similar  argument  we  find  the  following  relations 


where  Cn(z,u)  and  Dn(z,v)  are  defined  in  (2.19). 


(A. 7) 


Ol+,  = JZMtZ.] 
(2nq+p-q) 


+ qCp-qJzjfljj 


3nz2(2nq+i 


(2nq+p-q) 


(A. 8) 


^n  = C^(z,v)  or  D*(z,v)  , 


C*(z,v)  = 0 , c£(z,v)  = (p+q)  , 


D0(z>v)  = 1/P  > Di(z,u)  = (p+q)  + qz  , 


v„n,k  - N*p)tal1)vn,|,  * a^n^.q)p) 


(A. 9) 


5n(2nq+p+c 


(2nq+p-q)  Yn-l,k-2  ' 


(A. 10) 


Yn,k  “ cn,k  or  dn,k  * k °,1,  ...,n 
and  Yn,k  * 0 if  k < 0 or  k > n , 


A 

* * - 


c 


o,  o 


0 


> 


0 > cl,o=  (p+q) 
1/p  > 

(p+q)  , dlfl  = q 


> (A. 11) 


Tables  A.l  and  A. II  give  the  coefficients  of  the  polynomials  in  the 
rational  approximations  to  the  error  function,  see  (2.1),  (2.18),  (3.2)  and 
(3.3).  Table  A. Ill  gives  the  error  associated  with  the  approximation  (2.1) 
for  v = j . The  introduction  to  Table  A. Ill  shows  how  the  table  may  be  used 
to  get  Rn(z,v)  for  other  values  of  v and  Tn(z,v)  . Table  A. IV  gives 
coefficients  of  the  polynomials  in  the  rational  approximation  to  the  exponen- 
tial function,  see  (4.1).  Tables  A.V,  A. VI  and  A. VII  list  the  coefficients 
of  Un(z2)  , Wn(z2)  and  Vn(z2)  , respectively,  see  (4.7)  and  (4.8).  These 
coefficients  are  pertinent  to  the  evaluation  of  the  circular  functions. 

Most  of  the  tables  in  the  Appendices  were  typed  by  the  IBM  1620 
computer  directly  on  stencils,  while  the  balance  of  the  report  was  done  on 
an  ordinary  typewriter.  The  typewriters  have  different  type  sizes.  The  com- 
puter has  no  lower  case  characters,  etc.,  and  so  a slight  variance  in  notation  N 
is  introduced.  For  example,  in  Table  A.I,AN*(  Z , 1 /2  ) corresponds  to 
An(z,5)  , etc.  Also  in  Table  A.Vj(JN(  Z**2 ) corresponds  to  Un(z2)  , etc. 
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TABLE  A. I 


f 


TABLES  OF  THE  C (EFFICIENTS  OF  THE  POLYNOMIALS 
Ag(»,£)  AND  B*(z,^) 

Note:  For  the  definition  of  A*(z,£)  and  bJ(z,j)  , see  (2.1) 
and  (A.l).  The  sequence  of  numbers  given  is  for  the 
highest  power  to  the  lowest  power,  respectively. 

e.g.,  A$(z ,\)  = -128z3  + 1932z2  - 9240z  + 45045 

AN*(Zt 1/2)  BN*(Zt 1/2) 


N = 00 


1 


6 

15 


60 

420 

945 


280 

3780 

20790 

45045 


01 


N - 02 


- h 

15 


32 

-210 

945 


N = 03 


-128 
1537 
- 9240 
45045 
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TABLE  A. I (Continued) 


AN*(Z, 1/2) 


BN*(7# 1 / 2 ) 


N = 04 


2048 
-43560 
5 40540 
-22  52250 
114  86475 


5040 
1 10880 
10  81080 
54  05400 
114  86475 


N = 05 


-8192 
2 81424 
-38  91888 
453  33288 
1745  94420 
9166  20705 


22176 
7 20720 
108  10800 
<518  91800 
4364  86050 
9166  20705 


N *»  06 


65536 
-28  85792 
699  73904 
-7914  94704 
89625  13560 
-3  27946  51890 
17  56856  35125 


1 92192 

86  48640 
1837  83600 
23279  25600 
1 83324  14100 
8 43291  04860 
17  56856  35125 


N = 07 


-2  62144 
16 1 40480 
-4353  52320 
91454  22000 
-9  17867  80800 
102  39962  73300 
-361  41044  94000 
1965  16931  86125 


8 23680 
490  08960 
13967  55360 
2 44432  18800 
28  10970  16200 
210  82276  21500 
948  70242  96750 
1965  16931  86125 
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TABIE  A. I (Concluded) 
AN*(Z, 1/2) 


N = 08 


83  88608 
-6277  2595?- 
2 51360  15040 
-53  38931  08320 
1038  337°5  78000 
-9618  6809 6 04600 
1 06381  16578  08900 
-3  65521  4S326  19250 
20  10368  21294  05875 


2 

9 

20 


N - 09 

32467 
-14  34070 
473  05494 
- 8681  49968 
1 61017  72510 
-14  07934  64071 
154  87281  04783 
-521  20657  37253 
2892  69648  41756 


54432 

92960 

38720 

72000 

40000 

82400 

4 

26000 

45 

86000 

321 

37500 

1407 

23125 

2892 

N * 10 


2684  35456 
-3  05076  71040 

182  46049  11200 

-6185  20520  21760 

1 82600  62172  35200 

-30  32474  61076  56000  105 

544  83460  50675  52800  1 501 

-4560  40860  39327  81600  15010 

49985  79524  65547  67600  1 04137 

-J  65462  23889  48456  42750  4 51260 

9 25084  33563  93642  75375  9 25084 
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BN*(Z, 1/2) 


280  05120 
21283  89120 
7 82183  00160 
179  90209  03680 
2810  97016  20000 
30358  47774  96000 
20098  96368  46000 
74723  58203  18000 
10368  21294  05875 


1182  43840 
1 11740  47880 
51  40059  72480 
1499  18408  64000 
30358  47774  96000 
40197  92736  92000 
48711  91614  84000 
65891  40704  94000 
25774  90584  11250 
69648  41756  23125 


9932  48256 
11  42235  49400 
642  50746  50000 
23130  26876  16000 
58693  05698  25600 
16908  59875  61600 
07493  23289  72000 
74932  32897  20000 
07343  03224  32500 
65153  13972  07500 
33563  93642  75375 


1 


TABLE  A. II 

TABLES  OF  THE  COEFFICIENTS  OF  THE  POLYNOMIALS 
C*(z,£)  AND  D*(z,f) 


Note:  For  the  definition  of  C*(z,|)  and  D*(z,£)  , see  (2.18) 
and  (A. 7).  The  sequence  of  numbers  given  is  for  the 
highest  power  to  the  lowest  power,  respectively. 


e.g.,  C|(z,£)  = 84z2  - 42 Oz  + 3465 


N = 03 


CN*(Z, 1/2) 


DN*(Z,1/2) 


N = 01 


N = 02 


0 

84 

-420 

3465 


40 

420 

1890 

3465 


TABLE  A. II  (Continued) 


CN*(Zt 1/2) 


DN*(Z,1/2) 


N = 04 


0 

.744 
23100 
-?0090 
6 >5675 


560 
10080 
83160 
3 60360 
6 75675 


N = 05 


0 

5104 
-82368 
17  29728 
-61  26120 
436  48605 


2016 
55440 
7 20720 
54  05400 
229  72950 
436  48605 


N = 06 


0 

-25376 
15  53552 
-171  53136 
3026  30328 
.10184  67450 
70274  25405 


14784 
5 76576 
108  10800 
1225  22400 
8729  72100 
36664  82820 
70274  25405 


N - 07 


0 

1 58528 
-53  60576 
2061  87696 
-1^288  52640 
3 07868  16060 
-10  03917  91500 
67  76445  92625 


54912 
28  82880 
735  13440 
11639  62800 
1 22216  09400 
8 43291  04860 
35  13712  70250 
67  76445  92625 
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TABI£  A. II  (Concluded) 


CN*(Z,1/2) 


DN*(Z,1/2) 


N a 08 


0 

-33  70624 
3395  37600 
-73102  91040 
22  56425  02800 
- 197-  17857  23000 
2873  21307  27300 
-9170  79015  35250 
60920  24887  69875 


16  47360 
1120  20480 
37246  80960 
7 82183  00160 
112  43880  64800 
1124  38806  48000 
7589  61943  74000 
31442  70909  78000 
60920  24887  69875 


N « 09 


0 

62 

23360 

198  88896 

5320 

97280 

-11584  58880 

2 

23480 

85760 

7 

01345 

14560 

59 

96736 

34560 

-122 

38237 

44000 

1124 

38806 

48000 

3311 

97300 

72000 

15179 

23887 

48000 

-26551 

62101 

59200 

1 

46732 

64245 

64000 

3 

79059 

32634 

57000 

9 

74723 

98203 

18000 

-11 

91329 

31137 

22000 

40 

20736 

42588 

11750 

78 

18098 

60588 

00625 

78 

18098 

60588 

00625 

N ® 10 


0 

-1105  68960 
1 65206  97600 
-58  55525  91360 
2759  72254  27200 
-42544  70268  04800 
10  54145  93612  32800 
-81  01039  31733  09600 
1117  46689  40671  23600 
-3471  23578  10107  47750 
22563  03257  65698  60375 


472  97536 
49662  41280 
25  70029  86240 
856  67662  80000 
20238  98516  64000 
3 52158  34189  53600 
45  48711  91614  84000 
428  87855  20939  92000 
2814  51549  81168  22500 
11570  78593  67024  92500 
22563  03257  65698  60375 
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TABLE  A. Ill 

TABLES  OF  THE  ERRCR  FCR  THE  ERROR  FUNCTION 

Here  we  give  the  values  of  R*(z,£)  = 1 2z'2ezRn(z,5)  | , see  (2.15), 
for  n = 2(l)lO  , r = 1(1)10  and  0 = 0,  n/2,  tt  where  z = re*9  . 

It  is  pertinent  to  point  out  that  R^(z,u)  , see  (2.15),  and 
T*(z,v)  = I u_1zuezTn(z,u)  | , see  (2.22),  can  both  be  obtained  from  Rj(z,^) 
since 


Rn(*,v)  = 22-ur(v+l)ni-«n-lRn(z,i){l+0(l/n)}  , 


and 


Tn(z,u)  = •24'2ur(»H)n3/2'l,z”1nin(z,i){w(l/n)}  . 


Hence,  the  tables  given  here  essentially  include  the  values  of 
Rq(z,v)  and  T*(z, v)  for  admissible  u fixed,  but  otherwise  arbitrary  and 
the  values  of  n , r and  3 mentioned  above. 


n = 2 


1 

0.896  (-3) 

0.747  (-3) 

0.747 

(-3) 

2 

0.509  (-1) 

0.295  (-1) 

0.354 

(-1) 

3 

0.651  ( 

0.219 

0.377 

4 

0.477  ( 1) 

0.774 

0.230 

( 1) 

5 

0.268  ( 2) 

0.175  ( 1) 

0.108 

( 2) 

6 

0.132  ( 3) 

0.290  ( 1) 

0.443 

( 2) 

7 

0.609  ( 3) 

0.375  ( 1) 

0.171 

( 3) 

8 

0.275  ( 4) 

0.395  ( 1) 

0.642 

( 3) 

9 

0.125  ( 5) 

0.349  ( 1) 

0.243 

( 4) 

10 

0.578  ( 5) 

0.262  ( 1) 

0.938 

( 4) 

TABLE  A. Ill  (Continued) 


n = 3 


1 

0.521  ( 

- 5) 

0.456  ( 

- 5) 

0.456  ( 

- 5) 

2 

0.111  ( 

- 2) 

0.746  ( 

- 3) 

0.853  ( 

- 3) 

3 

0.294  ( 

- 1) 

0.132  ( 

- 1) 

0.197  ( 

- 1) 

4 

0.344 

0.905  ( 

- 1) 

0.202 

5 

0.264  ( 

1) 

0.358 

0.136  ( 

1) 

6 

0.160  ( 

2) 

0.973 

0.719  ( 

1) 

7 

0.838  ( 

2) 

0.200  ( 

1) 

0.330  ( 

2) 

8 

0.402  ( 

3) 

0.331  ( 

1) 

0.138  ( 

3) 

9 

0.183  ( 

4) 

0.454  ( 

1) 

0.552  ( 

3) 

10 

0.813  ( 

4) 

0.531  ( 

1) 

0.214  ( 

4) 

n = 4 


1 

0.178  ( 

- 7) 

0.160  ( 

- 7) 

0.160  ( 

- 7) 

2 

0.147  ( 

- 4) 

0.107  ( 

- 4) 

0.119  ( 

- 4) 

3 

0.833  ( 

- 3) 

0.443  ( 

- 3) 

0.608  ( 

- 3) 

4 

0.162  ( 

- 1) 

0.567  ( 

- 2) 

0.107  ( 

- 1) 

5 

0.181 

0.373  ( 

- 1) 

0.108 

6 

0.144  ( 

1) 

0.158 

0.765 

7 

0.923  ( 

1) 

0.484 

0.442  ( 

1) 

8 

0.513  ( 

2) 

0.116  ( 

1) 

0.221  ( 

2) 

9 

0.259  ( 

3) 

0.227  ( 

1) 

0.100  ( 

3) 

10 

0.123  ( 

4) 

0.375  ( 

1) 

0.428  ( 

3) 

n = 5 


1 

0.401  (-10) 

0.368  (-10) 

0.368 

<- 

10) 

2 

0.130  ( 

- 6) 

0.998  ( 

- 7) 

0.109 

(- 

6) 

3 

0.160  ( 

- 4) 

0.949  ( 

- 5) 

0.123 

(- 

4) 

4 

0.532  ( 

- 3) 

0.223  ( 

- 3) 

0.375 

(- 

3) 

5 

0.879  ( 

- 2) 

0.238  ( 

- 2) 

0.569 

(- 

2) 

6 

0.949  ( 

- 1) 

0.153  ( 

- 1) 

0.563 

(- 

1) 

7 

0.774 

0.678  ( 

- 1) 

0.421 

8 

0.520  ( 

1) 

0.227 

0.259 

( 

1) 

9 

0.306  ( 

2) 

0.609 

0.139 

( 

2) 

10 

0.161  ( 

3) 

0.135  ( 

1) 

0.677 

( 

2) 

30 


TABI£  A. Ill  (Continued) 

0 

nZf- 

rr 

0.639  (-13) 

0.593  (-13) 

0.593  (-13) 

0.812  (-  9) 

0.650  (-  9) 

0.700  (-  9) 

0.220  (-  6) 

0.141  (-  6) 

0.176  (-  6) 

0.127  (-  4) 

0.603  (-  5) 

0.941  (-  5) 

0.316  (-  3) 

0.104  (-  3) 

0.218  (-  3) 

0.470  (-  2) 

0.993  (-  3) 

0.302  (-  2) 

0.498  (-  1) 

0.625  (-  2) 

0.296  (-  1) 

0.414 

0.287  (-  1) 

0.229 

0.288  ( 1) 

0.103 

0.148  ( 1) 

0.176  ( 2) 

0.300 

0.841  ( 1) 

0.757  (-16) 

0.710  (-16) 

0.710  (-16) 

0.380  (-11) 

0.313  (-11) 

0.334  (-11) 

0.228  (-  8) 

0.155  (-  8) 

0.188  (-  8) 

0.228  (-  6) 

0.120  (-  6) 

0.176  (-  6) 

0.866  (-  5) 

0.329  (-  5) 

0.627  (-  5) 

0.180  (-  3) 

0.465  (-  4) 

0.122  (-  3) 

0.250  (-  2) 

0.411  (-  3) 

0.159  (-  2) 

0.261  (-  1) 

0.256  (-  2) 

0.156  (-  1) 

0.220 

0.121  (-  1) 

0.123 

0.158  ( 1) 

0.455  (-  1) 

0.829 

0.693  (-19) 

0.655  (-19) 

0.655  (-19) 

0.138  (-13) 

0.116  (-13) 

0.123  (-13) 

0.184  (-10) 

0.131  (-10) 

0.155  (-10) 

0.322  (-  8) 

0.182  (-  8) 

0.256  (-  8) 

0.187  (-6) 

0.793  (-  7) 

0.140  (-  6) 

0.547  (-  5) 

0.165  (-  5) 

0.388  (-  5) 

0.101  (-  3) 

0.203  (-  4) 

0.674  (-  4) 

0.133  (-  2) 

0.170  (-  3) 

0.840  (-  3) 

0.138  (-  1) 

0.104  (-  2) 

0.817  (-  2) 

0.117 

0.503  (-  2) 

0.659  (-  1) 

TABIE  A.  IV 


TABLE  OF  THE  COEFFICIENTS  OF  THE  POLYNOMIAL  Gn(z) 

For  the  definition  of  Gn(z)  , see  (4.1).  The  sequence 
of  numbers  given  is  for  the  highest  power  to  the 
lowest  power,  respectively. 

e.g.,  G4(z)  = z4  + 20z3  + 160z2  + 840z  + 1680 

GN  (Z ) 

1 


TABIE  A. IV  (Continued) 


30 

3360 

30240 


N = 06 


1 

840 
75600 
6 65280 


42 

10080 
3 32640 


N = 07 


N = 08 


N - 09 


1 

1512 
2 77200 
86  48640 


1 

2520 
8 31600 
605  40480 
5189  18400 


1 

3560 

21  62160 
3027  02400 
88216  12800 


- 34  - 


\ 


56 

25200 
19  °5840 
172  57280 


72 
55440 
86  48640 
2594  59200 


90 

1 10880 
302  70240 
20756  73600 
1 76432  25600 


TABIE  A. IV  (Concluded) 


r 


H * 

GN(  7.) 


N = 10 


1 110 
5940  2 05920 

50  45040  908  10720 


12108  09600  1 17621  50400 

7 93945  15200  33  52212  86400 

67  04425  72800 


TABLE  A.V 


; 1 


TABLE  OF  THE  COEFFICIENTS  OF  THE  POLYNOMIAL  Un(z2) 


Note:  For  the  definition  of  Uj^z2)  , see  (4.8).  The  sequence 
of  numbers  given  is  for  the  lowest  power  to  the 
highest  power,  respectively. 


e.g.,  U3(z  ) = -z6  + 264z4  - 6480z2  + 14400 


UN( Z**2 ) 


N = 00 

.10000  00000  00000  00000  +01 


M = 01 

.40000  00000  00000  00000  +01 


02 


. 14400  00000  00000  00000  +03 
. 10000  00000  00000  00000  +01 


03 


. 14400  00000  00000  00000  +05 
.26400  00000  00000  00000  +03 


N * 04 


.28224  00000  00000  00000  +07 
.6? 360  00000  00000  00000  +05 
.10000  00000  00000  00000  +01 


-.10000  00000  00000  00000  +01 


-.60000  00000  00000  00000  +02 


-.64800  00000  00000  00000  +04 
-.10000  00000  00000  00000  +01 


-.13104  00000  00000  00000  +07 
-.76000  00000  00000  00000  +03 
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1 


— 


TABI£  A.V  (Continued) 


UN (7**2) 


N = 05 


.51445 

76000 

00000 

00000 

+05 

-.43182 

72000 

00000 

00000 

+09 

.25804 

80000 

00000 

00000 

+08 

-.40824 

00000 

00000 

00000 

+06 

. 1 7400 

00000 

00000 

00000 

+04 

-. 10000 

00000 

00000 

00000 

+01 

I 

N 

= 06 

1 

• 

.44255 

74784 

00000 

00000 

+ 1? 

-.21123 

97056 

00000 

00000 

+ 12 

. 13535 

05280 

00000 

00000 

+ 1 1 

-.25788 

67200 

00000 

00000 

+0° 

. 17035 

20000 

00000 

00000 

+07 

-. 34440 

00000 

00000 

00000 

+04 

. 10000 

00000 

00000 

00000 

+01 

1 

'I 

N 

= 07 

. 299  1c 

58553 

98400 

00000 

+ 15 

-. 14384 

41804 

80000 

00000 

+15 

. 5641'5 

66233 

60000 

00000 

+ 13 

-.20552 

09587 

20000 

00000 

+ 12 

N 

. 17141 

24160 

00000 

00000 

+ 10 

-. 56629 

44000 

00000 

oooco 

+07 

1 

! 

.61600 

00000 

00000 

00000 

+04 

-. 10000 

00000 

00000 

00000 

+01 

= 08 

• 

. 26527 

63058 

58560 

00000 

+ 18 

-.13015 

02144 

98304 

00000 

+ 18 

a 

.50161 

53232 

48640 

00000 

+ 16 

-. 20687 

40850 

17600 

00000 

+15 

i ! 

i 

. 15540 

43744 

00000 

00000 

+ 13 

-.86313 

42720 

00000 

00000 

+ 10 

. 15556 

56000 

00000 

00000 

+08 

-. 10224 

00000 

00000 

00000 

+05 

. 10000 

00000 

00000 

00000 

+01 

N 

=>  05 

.31128 

34055 

72495 

36000 

+21 

-. 15106 

40075 

86652 

16000 

+21 

. 10717 

15657 

31706 

88000 

+20 

-.  25935 

10574 

05440 

00000 

+ 18 

.27586 

13803 

54560 

00000 

+ 16 

-.14176 

33176 

96000 

00000 

+14 

.35472 

72560 

00000 

00000 

+ 11 

-. 39564 

32000 

00000 

00000 

+08 

r | 

| 

F I 

. 16020 

00000 

00000 

00000 

+05 

- 37  - 

-. 10000 

00000 

00000 

00000 

+01 

Ll 

* . v 

\ 

TABI£  A.V  (Concluded) 


UN (7**2) 


N = 10 

b44c*1»q  37434  7.7.683  29984  +24 
.15812  89202  64694  57520  +23 
.4545^  37982  88175  20000  +19 
.79582  79565  28000  00000  +14 
.90676  08000  00000  00000  +08 
.10000  00000  00000  00000  +01 


-.21883  22369  29^64  23808  +24 
-.39825  96563  64810  24000  +21 
-.26326  17302  51520  00000  +17 
-.12473  80992  00000  00000  +12 
-.23980  00000  00000  00000  +05 


f 


I 


K 

I 


TABLE  A. VI 

TABIES  OF  THE  COEFFICIENTS  OF  THE  POLYNOMIAL  Vin(z2) 

Note:  For  the  definition  of  Wn(z2)  , see  (4.8).  The  sequence 
of  numbers  given  is  for  the  lowest  power  to  the 
highest  power,  respectively. 

e.g.,  W3(z2)  = z6  + 24z4  + 720z2  + 14400 

V/N(Z**2 ) 

N - 00 

. 10000  00000  00000  00000  +01 

N * 01 

.40000  00000  00000  00000  +01  . 10000  00000  00000  00000  +01 


N a 0? 

14400  00000  00000  00000  +03  .12000  00000  00000  00000  +02 

! 10000  00000  00000  00000  +01 


N - 03 


. 14400  00000  00000 
.24000  00000  00000 


00000  +05 
00000  +02 


.72000  00000  00000  00000  +03 
.10000  00000  00000  00000  +01 


N - 04 


' 





.28224  00000  00000  00000  +07 
.21600  00000  00000  00000  +04 
. 10000  00000  00000  00000  +01 


.10080  00000  00000  00000  +06 
.40000  00000  00000  00000  +02 


A 


- 39 


TABUS  A. VI  (Continued) 


WN(Z**2 ) 


. 9 1445  76000  00000  00000  +09 
.40320  00000  00000  00000  +06 
.60000  00000  00000  00000  +02 


. 4425S  74784  00000  00000  +12 
.12700  80000  00000  00000  +09 
. 10080  00000  00000  00000  +05 
.10000  00000  00000  00000  +01 


. 299 11*  58953  ^8400  00000  +15 
.60354  20160  00000  00000  +11 
.30240  00000  00000  00000  +07 
.11200  00000  00000  00000  +03 


. 25401  60000  00000  00000  +08 
.50400  00000  00000  00000  +04 
.10000  00000  00000  00000  +01 


.10059  03360  00000  00000  +11 
.12096  00000  00000  00000  +07 
.84000  00000  00000  00000  +02 


.57537  67219  20000  00000  +13 
.46569  60000  00000  00000  +09 
.18144  00000  00000  00000  +05 
. 10000  00000  00000  00000  +01 


N = 08 

.26927  63058  58560  00000  +18 
.40276  37053  44000  00000  +14 
. 13t'70  88000  00000  00000  +10 
. 30240  00000  00000  00000  +05 
. 10000  00000  00000  00000  +01 


.44879  38430  97600  00000  +16 
.26153  48736  00000  00000  +12 
.66528  00000  00000  00000  +07 
. 14400  00000  00000  00000  +03 


.31128  34095  72455  36000  +21 
.35903  50744  78080  00000  +17 
.91537  20576  00000  00000  +12 
.13305  60000  00000  00000  +08 
. 18000  00000  00000  00000  +03 


.45776  97199 
.20138  18526 
.36324  28800 
.47520  00000 
. 10000  00000 


59552  00000  +19 
72000  00000  +15 
00000  00000  +10 
00000  00000  +05 
00000  00000  +01 


TABLE  A. VI  (Concluded) 


1 


WN(Z**2) 


N = 10 

J'MM  32434  2? 68 3 
.41  If S 2747?  63556 
.80552  74106  88000 
.84756  67200  00000 
.71280  00000  00000 
. 10000  00000  00000 


2C!?84  +24  .55143 
80000  +20  .20345 
00000+15  .27461 
00000  +10  .24710 
00000  +05  .22000 
00000  +01 


84781  87741  18400  +22 
32088  70S  12  00000  +18 
16172  80000  00000  +13 
40000  00000  00000  +08 
00000  00000  00000  +03 
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TABlE  A. VII 


TABLES  OF  THE  COEFFICIENTS  OF  THE  POLYNOMIAL  VQ(z2) 


Note:  For  the  definition  of  Vn(z2)  , see  (4.8).  The  sequence 
of  numbers  given  is  for  the  lowest  power  to  the 
highest  power,  respectively. 

e.g.,  V3(z2)  = 24 Z-  - 1680z2  + 14400 


VN(Z**2) 


N - 00 

.00000  00000  00000  00000  00 


N = 01 

.40000  00000  00000  00000  01 


N = 02 

.14400  00000  00000  00000  03  -.12000  00000  00000  00000  (U 


N - 03 

. 14400  00000 
. 7.4000  00000 

N = 04 

. 28224  00000 
.88800  00000 


00000  00000  +05 
0C000  00000  +02 


00000  00000  +07 
00000  00000  +04 


-.  16800  00000 


-.36560  00000 
-.40000  00000 


00000  00000  +04 


00000  00000  +06 
00000  00000  +0^ 
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TABLE  A. VII  (Continued) 


N = 05 

.91445  76000 
. 37900  80000 
.60000  00000 


N - 06 

.44259  74784 
.21388  14720 
.90720  00000 


N = 07 

.29919  58953 
. 15946  92126 
.11124  28800 
.11200  00000 


N - 08 

.26927  63058 
.15362  55847 
.14503  37011 
.47376  00000 


N - 09 

.31128  34095 
.18669  82387 
.21608  80001 
. 13278  98880 
. 18000  00000 


00000  00000  +09 
00000  00000  +07 
00000  00000  +02 


00000  00000  +12 
00000  00000  +10 
00000  00000  +05 


98400  00000  +15 
77000  00000  +13 
00000  00000  +09 
00000  00000  +03 


58560  00000  +18 
52640  00000  +16 
20000  00000  +12 
00000  00000  +06 


72495  36000  +21 
28601  60000  +19 
28000  00000  +15 
00000  00000  +10 
00000  00000  +03 


-.12700  80000 
-.31920  00000 


-.63707  21280 
-.23950  08000 
-.84000  00000 


-.44112  21534 
-.21009  54240 
-.21974  40000 


-.40391  44587 
-.22479  54508 
-.41646  52800 
-. 14400  00000 


-.47302  87106 
-.29397  64065 
-.77785  86816 
-.93456  00000 


00000  00000  +09 
00000  00000  +05 


00000  00000  +11 
00000  00000  +08 
00000  00000  +02 


72000  00000  +14 
00000  00000  +11 
00000  00000  +06 


87840  00000  +17 
80000  00000  +14 
00000  00000  +09 
00000  00000  +03 


24870  40000  +20 
63840  00000  +17 
00000  00000  +12 
00000  00000  +06 
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TABUS  A. VII  (Concluded) 


N * 10 

. 44s 49  32434  22683 
.28012  45506  33915 
.37541  77850  14272 
.34464  83040  00000 
.17186  40000  00000 


29984  +24  -.6900 1 
18720  +22  -.46561 
00000  +18  -.15728 
00000  +13  -.37378 
00000  +07  -.22000 


15578  85698  04800  +23 
72008  73144  32000  +20 
18837  99040  00000  +16 
59840  00000  00000  +10 
00000  00000  00000  +03 
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APPENDIX  B 
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COEFFICIENTS  OF  RATIONAL  APPROXIMATIONS  TO  E(z) 
AND  RELATED  INTEGRALS 


Table  B.I  gives  the  coefficients  of  the  polynomials  of  the  main 
diagonal  Pade'  approximants  of  E(z),  see  (5.1).  Displayed  in  Table  B.II  is 
the  absolute  value  of  the  errors  associated  with  the  approximations  listed 
in  B.I.  Similar  coefficients  for  the  main  diagonal  Fade'  approximants  to 
z-1Si(z)  , see  (5.2),  are  given  for  n even  in  Table  B.III.  Table  B.IV 
gives  the  error  associated  with  the  approximations  listed  in  B.III.  Table 
B.V  lists  the  coefficients  of  the  polynomials  of  the  main  diagonal  Pade' 
approximants  to  4z"2H(z)  , see  (5.3)  , for  even  n . Table  B.VI  lists  the 
errors  associated  with  the  approximants  in  B.V.  Tables  B.VII  and  B.VIII 
Give  the  coefficients  which  define  the  approximations  (5.11)  to  E(z)  , and 
the  corresponding  error  tables,  respectively. 


TABLE  B.I 


COEFFICIENTS  OF  PAPE  APPROXIMATION  TO  E(z) 

Note:  For  the  definition  of  E(z)  and  its  Pade'  approximants,  see 
Section  V.  The  coefficients  are  given  for  n = l(l)lO  . The 
expression  attached  to  the  numbers  on  the  right  indicates 
the  power  of  10  by  which  the  number  is  multiplied. 

e.g.,  E2(z)  = ^ + 0 • Q96z  + 0 . 04555 . . . z^ 

1 + 0.336Z  + 0.0333. . . z2 


NUMERATOR 


DENOMINATOR 


. 10000  00000  00000  00000  +01 
.27777  77777  77777  77778  -01 


M - 02 

. 10000  00000  00000  00000  +01 

.86000  00000  00000  00000  -01 

.45555  55555  55555  55556  -02 


. 10000  00000  00000  00000  +01 

.11548  60120  33^92  92848  +00 

. 14987  35259  04789  47405  -01 

-.85008  16740  06988  81376  -04 


. 10000  00000  00000  00000  +01 
.15183  03075  00550  75001  +00 
.22174  10124  06340  46892  -01 
.71679  28213  71348  661  18  -03 
. 10529  22995  44767  84505  -04 


.10000  00000  00000  00000  +01 
.22222  22222  22222  22222  +00 


. 10000  00000  00000  00000  +01 
.33600  00000  00000  00000  +00 
.33000  00000  00000  00000  -01 


. 10000  00000  00000  00000  +01 
.36548  60120  33992  92848  +00 
.50803  30004  34216  23969  -01 
.27277  05063  78843  33065  -02 


. 10000  00000  00000  00000  +01 
.40183  03075  00550  75001  +00 
.67076  12256  02161  78840  -01 
.55785  84155  83924  05926  -02 
. 19778  97186  99680  64309  -03 


COLUMBIA  UNIVERSITY 
hi!os;n  LABORATORIES 
CONTRACT  Nonr-266(84) 


r 


f 


TABIE  B. I (Continued) 


NUMERATOR 


DENOMINATOR 


N = 05 


. 10000  00000 
* 16  7. 04  16009 
. 265? 5 51845 
. 12747  8327?- 
.60877  80145 
-.  13654  42246 


00000  00000  +01 
62342  70475  +00 
2?<->?3  74506  -01 
48448  67889  -02 
04733  48711  -04 
45193  12376  -06 


. 10000  00000  00000  00000  +01 
.41204  16009  62342  70975  +00 
.74050  36313  73224  96387  -01 
.73128  40670  01829  36777  -02 
.40061  22491  28200  19535  -03 
.98510  21164  96373  17832  -05 


N = 06 


. 10000  00000  00000 
. 17964  76037  82301 
. 30906  3‘'014  59549 
. 19643  87594  46735 
. 11565  07429  44841 
.19803  30657  69 150 
.12947  26971  80815 


00000  +01 
29344  +00 
31654  -01 
53361  -02 
78517  -03 
50495  -05 
10906  -07 


. 10000  00000 
.42964  76037 
.82762  73553 
.92024  26823 
.62715  69030 
.25037  59209 
.46180  79516 


00000  00000  +01 
82301  29344  +00 
59746  99458  -01 
88873  61206  -02 
34069  87432  -03 
97727  45364  -04 
66610  03704  -06 


N - 07 


. 10000  00000 
.18480  11227 
.33261  21378 
.23552  22731 
.16377  94153 
.43980  12290 
. 11014  77035 
12839  85657 


00000  00000  +01 
26475  81023  +00 
35356  62215  -01 
54871  76423  -02 
15153  96938  -03 
60812  99224  -05 
81094  26833  -06 
45606  74048  -09 


. 10000  00000  00000 
.43480  11227  26475 
.86405  93890  95990 
. 10217  75619  63664 
.78040  02200  52063 
.38719  75252  89419 
. 1 1633  17088  02093 
.16457  10741  45580 


00000  +01 
81023  +00 
59216  -01 
92989  -01 
61676  -03 
77095  -04 
24466  -05 
73556  -07 
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TABLE  B.I  (Concluded) 


NUMERATOR 


DENOMINATOR 


N = 08 

. 10000 
. 19513 
. 36078 
. 28530 
.71463 
.77757 
. 25482 
. 27089 
.97825 


N = 09 

. 10000 
. 19823 
37482 


00000 

36973 

77455 

41473 

38297 

97735 

33856 

12491 

79525 


00000 

48381 

48033 

63713 

93936 

47219 

33725 

61763 

95621 


00000 

21416 

66895 

87721 

23562 

15163 

45105 

77927 

99712 


+01 

+00 

-01 

-02 

-03 

-05 

-06 

-08 

-11 


. 10000 
.44513 
.91756 
.11479 
. 95698 
. 54679 
.21015 
. 49960 

. 56625 


00000 

36973 

59333 

20657 

96347 

34064 

17140 

40804 

36991 


00000 

48381 

63431 

17017 

45946 

85047 

91386 

95845 

72465 


00000 

21416 

14880 

38378 

53974 

10025 

35547 

41259 

57055 


+01 

+00 

-01 

-01 

-03 

-04 

-05 

-07 

-09 


00000 

84189 

24827 


00000 

19627 

32679 


00000 

19265 

61682 


+01 

+00 

-01 


10000 

44823 


00000 

84189 

29744 


00000 
19627 
76 1 S 2 


00000  +01 
19265  +00 
04289  -01 


H 

.31215 

54261 

10534 

67154  -02 

. 12132 

66090 

52530 

81529  -01 

! 

.25168 

68581 

42141 

13445  -03 

. 10658 

74645 

55601 

67575  -02 

. 

. 10438 

66937 

291 14 

63287  -04 

.66311 

99549 

82139 

94997  -04 

.41779 

8<?027 

92798 

76366  -06 

. 29328 

17107 

87334 

24422  -05 

. 72282 

10040 

34696 

94233  -08 

.89307 

30009 

59419 

32977  -07 

j 

. 11156 

38762 

07c-00 

30789  -09 

. 17075 

36697 

26887 

18576  -08 

- . 78469 

* 

72732 

66716 

49614  -13 

. 15708 

67540 

09828 

56485  -10 

n 

..  • 

N =*  10 

• 

; 

. 10000 

00000 

00000 

00000  +01 

. 10000 

00000 

00000 

00000  +01 

■ 

• 

. 20502 

08456 

77917 

06628  +00 

. 45502 

08456 

77917 

06628  +00 

. 39390 

75193 

16296 

55245  -01 

.97590 

40779 

55533 

66258  -01 

. 34858 

23655 

29  237 

91179  -02 

. 13021 

15639 

98519 

94778  -01 

76100  -02 

. 29317 

75506 

142.66 

48946  -03 

. 1 1999 

11137 

74704 

I 

. 13754 

73570 

2.9922 

39386  -04 

.80015 

09559 

21661 

45984  -04 

; 

. 60^64 

46174 

77455 

80030  -06 

.39222 

83073 

88575 

92254  -05 

■ 

. 14447 

18655 

00891 

74805  -07 

. 14003 

621 18 

96032 

45150  -06 

. 30430 

04327 

31332 

24684  -09 

. 34984 

41348 

05290 

45646  -08 

. 22059 

38908 

74765 

26250  -11 

. 55465 

89453 

73869 

45817  -10 

. 49848 

28058 

16872 

88340  -14 

.42591 

33901 

24021 

43020  -12 
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] 

TABUS  B. II 


ERROR  OF  PAPE  APPROXIMATIONS  TO  E(z) 


Let  En(z)  be  the  n^h  order  main  diagonal  Pade*  approximation  to 
E(z),  see  B.I,  and  define  en(z)  = |E(z)  - En(z)  | . The  tables  give  en(z) 
for  n = 2(1)10  , r = l(l)lO  and  9 = 0,  rr/2  , n where  z = re^®  . The 
expression  in  parentheses  to  the  right  of  each  number  indicates  the  power  of 
10  by  which  the  number  is  multiplied. 


* 

rx 


n = 2 


tt/2 


rr  (Relative  Error) 


1 

0.782  (-  5) 

0.107  (-  4) 

0.320  ( 

- 4) 

0.243  ( 

- 4) 

2 

0.137  (-  3) 

0.418  (-  3) 

0.232  ( 

- 2) 

0.126  ( 

- 2) 

3 

0.599  (-  3) 

0.267  (-  2) 

0.433  ( 

- 1) 

0.157  ( 

- 1) 

4 

0.152  (-  2) 

0.849  (-  2) 

0.455 

0.103 

5 

0.292  (-  2) 

0.202  (-  1) 

0.288  ( 

1) 

0.379 

6 

0.473  (-  2) 

0.361  (-  1) 

0.102  ( 

1) 

0.732 

7 

0.686  (-  2) 

0.542  (-  1) 

0.247  ( 

2) 

0.915 

8 

0.922  (-  2) 

0.713  (-  1) 

0.533  ( 

2) 

0.974 

9 

0.117  (-  1) 

0.845  (-  1) 

0.114  ( 

3) 

0.991 

10 

0.143  (-  1) 

0.825 

0.248  ( 

3) 

0.996 

n = 3 

1 

0.217  (-  7) 

0.443  (-  7) 

0.101  ( 

- 6) 

0.776  ( 

- ?) 

2 

0.139  (-  5) 

0.523  (-  5) 

0.302  ( 

- 4) 

0.164  ( 

- 4) 

3 

0.125  (-  4) 

0.780  (-  4) 

0.128  ( 

- 2) 

0.465  ( 

- 3) 

4 

0.510  (-  4) 

0.482  (-  3) 

0.252  ( 

- 1) 

0.571  ( 

- 2) 

5 

0.138  (-  3) 

0.179  (-  2) 

0.345 

0.454  ( 

- 1) 

6 

0.291  (-  3) 

0.474  (-  2) 

0.454  ( 

1) 

0.326 

7 

0.519  (-  3) 

0.985  (-  2) 

0.233  ( 

3) 

0.863 

6 

0.823  (-  3) 

0.170  (-  1) 

0.703  ( 

2) 

0.128  ( 

1) 

9 

0.120  (-  2) 

0.255  (-  1) 

0.123  ( 

3) 

0.107  ( 

1) 

10 

0.164  (-  2) 

0.340  (-  1) 

0.254  ( 

3) 

0.102  ( 

1) 
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TABLE  B. II  (Continued) 


V 

rj 

0 

n/2 

n 

n (Relative  Error) 

n = 4 

1 

0.593  (-10) 

0.128  (-  9) 

0.302  (-  9) 

0.229  (-  9) 

2 

0.144  (-  7) 

0.610  (-  7) 

0.375  (-  6) 

0.204  (-  6) 

3 

0.273  (-  6) 

0.209  (-  5) 

0.366  (-  4) 

0.133  (-  4) 

4 

0.186  (-  5) 

0.236  (-  4) 

0.130  (-  2) 

0.294  (-  3) 

5 

0.738  (-  5) 

0.141  (-  3) 

0.271  (-  1) 

0.366  (-  2) 

6 

0.209  (-  4) 

0.558  (-  3) 

0.401 

0.288  (-  1) 

7 

0.474  (-  4) 

0.162  (-  2) 

0.428  ( 1) 

0.159 

8 

0.920  (-  4) 

0.375  (-  2) 

0.273  ( 2) 

0.499 

9 

0.159  (-  3) 

0.720  (-  2) 

0.945  ( 2) 

0.822 

10 

0.251  (-  3) 

0 119  (-  1) 

0.237  ( 3) 

0.952 

n = 5 

1 

0.837  (-13) 

0.187  (-12) 

0.451  (-12) 

0.342  (-12) 

2 

0.781  (-10) 

0.362  (-  9) 

0.227  (-  8) 

0.123  (-  8) 

3 

0.319  (-  8) 

0.284  (-  7) 

0.502  (-  6) 

0.182  (-  6) 

4 

0.368  (-  7) 

0.588  (-  6) 

0.316  (-  4) 

0.715  (-  5) 

5 

0.206  (-  6) 

0.574  (-  5) 

0.102  (-  2) 

0.134  (-  3) 

6 

0.834  (-  6) 

0.342  (-  4) 

0.217  (-  1) 

0.156  (-  2) 

7 

0.243  (-  5) 

0.143  (-  3) 

0.355 

0.131  (-  1) 

8 

0.579  (-  5) 

0.457  (-  3) 

0.513  ( 1) 

0.937  (-  1) 

9 

0.119  (-  4) 

0.117  (-  2) 

0.104  ( 3) 

0.904 

10 

0.218  (-  4) 

0.250  (-  2) 

0.442  ( 3) 

0.178  ( 1) 

n = 6 

1 

0.117  (-15) 

0.268  (-15) 

0.653  (-15) 

0.499  (-15) 

2 

0.423  (-12) 

0.209  (-11) 

0.135  (-10) 

0.733  (-11) 

3 

0.376  (-10) 

0.373  (-  9) 

0.677  (-  8) 

0.246  (-  6) 

4 

0.746  (-  9) 

0.139  (-  7) 

0.763  (-  6) 

0.173  (-  6) 

5 

0.658  (-  8) 

0.217  (-  6) 

0.384  (-  4) 

0.505  (-  5) 

6 

0.352  (-  7) 

0.192  (-  5) 

0.118  (-  2) 

0.847  (-  4) 

7 

0.134  (-  6) 

0.113  (-  4) 

0.258  (-  1) 

0.956  (-  3) 

8 

0.399  (-  6) 

0.489  (-  4) 

0.445 

0.813  (-  2) 

9 

0.993  (-  6) 

0.165  (-  3) 

0.620  ( 1) 

0.539  (-  1) 

10 

0.215  (-  5) 

0.455  (-  3) 

0.627  ( 2) 

0.252 
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TABLE  B. II  (Continued) 

r' 

0 

n/2  tt 

n (Relative  Error) 

n = 7 


1 

0.900  (-19) 

0.233  (-18) 

0.580  (-18) 

0.440  (-18) 

2 

0.140  (-14) 

0.722  (-14) 

0.472  (-13) 

0.256  (-13) 

3 

0.272  (-12) 

0.294  (-11) 

0.536  (-10) 

0.195  (-10) 

4 

0.932  (-11) 

0.198  (-  9) 

0.107  (-  7) 

0.242  (-  8) 

5 

0.125  (-  9) 

0.493  (-  8) 

0.840  (-  6) 

0.110  (-  6) 

6 

0.925  (-  9) 

0.644  (-  7) 

0.367  (-  4) 

0.263  (-  4) 

7 

0.462  (-  8) 

0.534  (-  6) 

0.108  (-  2) 

0.400  (-  3) 

8 

0.173  (-  7) 

0.314  (-  5) 

0.241  (-  1) 

0.440  (-  3) 

9 

0.524  (-  7) 

0.141  (-  4) 

0.436 

0.379  (-  2) 

10 

0.134  (-  6) 

0.336 

0.686  ( 1) 

0.276  (-  1) 

n = 8 

1 

0.370  (-20) 

0.100  (-19) 

0.100  (-19) 

0.759  (-20) 

2 

0.459  (-17) 

0.246  (-16) 

0.164  (-15) 

0.890  (-16) 

3 

0.197  (-14) 

0.227  (-13) 

0.421  (-12) 

0.153  (-12) 

4 

0.118  (-12) 

0.275  (-11) 

0.150  (-  9) 

0.340  (-10) 

5 

0.239  (-11) 

0.108  (-  9) 

0.184  (-  7) 

0.242  (-  8) 

6 

0.249  (-10) 

0.207  (-  8) 

0.115  (-  5) 

0.825  (-  7) 

7 

0.165  (-  9) 

0.238  (-  7) 

0.461  (-  4) 

0.171  (-  5) 

8 

0.784  (-  9) 

0.187  (-  6) 

0.133  (-  2) 

0.243  (-  4) 

9 

0.291  (-  8) 

0.109  (-  5) 

0.301  (-  1) 

0.262  (-  3) 

10 

0.895  (-  8) 

0.500  (-  5) 

0.565 

0.227  (-  2) 

n = 9 

1 

0.367  (-20) 

0.100  (-19) 

0.100  (-19) 

0.758  (-20) 

2 

0.500  (-20) 

0.500  (-19) 

0.377  (-18) 

0.205  (-18) 

3 

0.950  (-17) 

0.117  (-15) 

0.217  (-14) 

0.788  (-15) 

4 

0.998  (-15) 

0.254  (-13) 

0.138  (-11) 

0.312  (-12) 

5 

0.311  (-13) 

0.106  (-10) 

0.263  (-  9) 

0.346  (-10) 

6 

0.456  (-12) 

0.442  (-10) 

0.236  (-  7) 

0.169  (-  8) 

7 

0.400  (-11) 

0.707  (-  9) 

0.127  (-  5) 

0.470  (-  7) 

8 

0.242  (-10) 

0.743  (-  8) 

0.474  (-  4) 

0.866  (-  6) 

9 

0.111  (-  9) 

0.564  (-  7) 

0.134  (-  2) 

0.117  (-  4) 

10 

0.408  (-  9) 

0.329  (-  6) 

0.306  (-  1) 

0.123  (-  3) 
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TABIE  B.II  (Concluded) 


c 

> 

nZi 

T 

r 

tt  (Relative  Error) 

0.367 

(-20) 

0.100 

(-19) 

0.100 

(-19) 

0.759 

(-20) 

0.551 

(-20) 

0.141 

(-19) 

0.200 

(-20) 

0.109 

(-20) 

0.145 

(-18) 

0.578 

(-18) 

0.111 

(-16) 

0.403 

(-17) 

0.850 

(-17) 

0.231 

(-15) 

0.126 

(-13) 

0.285 

(-14) 

0.407 

(-15) 

0.227 

(-13) 

0.375 

(-11) 

0.493 

(-12) 

0.045 

(-14) 

0.923 

(-12) 

0.485 

(-  9) 

0.341 

(-10) 

0.989 

(-13) 

0.203 

(-10) 

0.354 

(-  7) 

0.131 

(-  8) 

0.766 

(-12) 

0.284 

(-  9) 

0.172 

(-  5) 

0.314 

(-  7) 

0.434 

(-11) 

0.277 

(-8) 

0.610 

(-  4) 

0.530 

(-  6) 

0.193 

(-10) 

0.204 

(-  7) 

0.170 

(-  2) 

0.683 

(-  5) 

TABIE  B.III 


COEFFICIENTS  OF  PAPE  APPROXIMATIONS  TO  z^SiCz) 


Note:  For  the  definition  of  Si(z)  and  its  Pade'  approximants,  see 
Section  V.  The  coefficients  are  given  for  n = 2(2)10.  The 
expression  attached  to  the  numbers  on  the  right  indicates 
the  power  of  10  by  which  the  numbers  are  multiplied. 


_1qj  ( _\  _ 1 “ 0.02555. . . z^ 
e.g.,  z Sio(z)  - ~ 


NUMERATOR 


N = 02 


.10000  00000  00000  00000  +01 
-.25555  55555  55555  55556  -01 


N = 04 


. 10000  00000  00000  00000  +01 
-.30427  89735  63441  27723  -01 
.51431  13199  43054  89416  -03 


06 


.10000  00000  00000  00000  +01 
-.35595  17688  83526  52167  -01 
.74030  23630  09289  15834  -03 
-.44148  31037  35430  42004  -05 


N 


08 


. 10000  00000  00000  00000  +01 
-.39097  46707  47744  00193  -01 
.88244  03534  29201  16961  -03 
“.75392  62090  91062  85582  -05 
.24448  76200  37211  27450  -07 


1 + 0. 03z 


DENOMINATOR 


.10000  00000  00000  00000  +01 

.30000  00000  00000  00000  -01 


.10000  00000  00000  00000  +01 

.25127  65819  92114  27833  -01 

.24362  56643  43689  77376  -03 


.10000  00000  00000  00000  +01 
.19960  37866  72029  03389  -01 
.18254  56222  98339  34660  -03 
.80396  58468  07101  69744  -06 


.10000  00000  00000  00000  +01 
.16458  08848  07811  55363  -01 
.13011  19356  94820  91197  -03 
.60370  25148  69333  23738  -06 
. 14412  35225  57797  97185  -08 
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TABI£  B.III  (Concluded) 


N = 10 


. 10000  00000  00000 
-.41604  56609  13051 
.98695  05591  58116 
-.99826  03983  42126 
.48154  24911  89950 
-.89943  23724  44804 


00000  +01 

. 10000 

18594  -01 

.13950 

83815  -03 

.95338 

42340  -05 

. 40702 

63472  -07 

.11122 

92329  -10 

. 16036 

00000  00000  00000  +01 
98946  42504  36962  -01 
86272  75855  58234  -04 
15961  74597  46660  -06 
28953  69313  70588  -08 
02483  25839  25315  -11 


Let  z"^-Sin(z)  be  the  nth  order  main  diagonal  (see  B.IIl)  Pade* 
approximant  to  z'^-Si(z)  and  define  en(z)  = | z"^Si(z  )-z-1Sin(z  ) |.  The 
tables  give  en(z ) for  n = 2(2)10  and  z = l(l)lO  , z real.  The  expres- 
sion in  parentheses  to  the  right  of  each  number  indicates  the  power  of  10  by 
which  the  number  is  multiplied. 


5_ 

2 

4 

e 

8 

10 

1 

0.210 

(-4) 

0.142 

(-8) 

0.251 

(-13) 

0.170 

(- 

18) 

0.100 

(-20) 

2 

0.112 

(-2) 

0.127 

(-5) 

0.369 

(- 

9) 

0.392 

(- 

13) 

0.186 

(-17) 

3 

0.990 

(-2) 

0.590 

(-4) 

0.898 

(- 

7) 

0.498 

(- 

10) 

0.122 

(-13) 

4 

0.401 

(-D 

0.775 

(-3) 

0.392 

(- 

5) 

0.714 

(- 

8) 

0.570 

(-11) 

5 

0.104 

0.496 

(-2) 

0.647 

(- 

4) 

0.302 

(- 

6) 

0.616 

(- 

9) 

6 

0.234 

0.198 

(-1) 

0.565 

(- 

3) 

0.579 

(- 

5) 

0.254 

(- 

7) 

7 

0.310 

0.563 

(-D 

0.313 

(- 

2) 

0.630 

(- 

4) 

0.539 

(- 

6) 

Q 

0.415 

0.125 

0.123 

(- 

1) 

0.449 

(- 

3) 

0.692 

(- 

5) 

9 

0.497 

0.227 

0.370 

(- 

1) 

0.229 

(- 

2) 

0.599 

(- 

4) 

10 

0.555 

0.335 

0.896 

(- 

1) 

0.889 

(- 

2) 

0.376 

(- 

3) 

TABLE  B.  V 


COEFFICIENTS  OF  PAPE  APPROXIMATIONS  TO  4z"2H(z' 


Note:  For  the  definition  of  H(z)  and  its  Pade7  approximants,  see 
Section  V.  The  coefficients  are  given  for  n = 2(2)10.  The 
expression  attached  to  the  numbers  on  the  right  indicates 
the  power  of  10  by  which  the  numbers  are  multiplied. 


e.g. , 4z"‘;H2(z)  = 


_ 1 - 0. 019444... z* 


1 + 0. 0222. . .z2 


NUMERATOR 


.10000  00000  00000  00000  +01 
1?444  44444  44444  44444  -01 


N = 04 


.10000  00000  00000  00000  +01 
■.21441  514756  55430  71161  -01 
.23504  84513  40586  17205  -03 


. 10000  00000  00000  00000  +01 

-.2466?  60143  3836?  60391  -01 

.34682  10542  75433  98340  -03 

15875  89095  57136  15182  -05 


. 10000  00000  00000  00000  +01 
-.27151  18247  38415  96633  -01 
.4215  7 85)824  50745  0?333  -03 
-.27389  37015  8??,04  22581  -05 
.70221  55017  2?417  74512  -08 


DENOMINATOR 


.10000  00000  00000  00000  +01 

.22222  22222  22222  22222  -01 


. 10000  00000  00000  00000  +01 
.20224  71910  11235  95506  -01 
.15181  c)  1 546  28143  39?- 19  -03 


.10000  00000  00000  00000  +01 
.17004  06523  28304  06275  -01 
.12930  78463  84108  31895  -03 
. 460?7  66426  80130  08786  -06 


.10000  00000  00000  00000  +01 


. 14475  48419  28250  700*4  -01 
.99198  23122  58636  82145  -04 


.39189  0164?  14937  8528?  -06 
.77891  33652  63414  35304  -09 


-J 


TABLE  B.V  (Concluded) 


NUMERATOR 


DENOMINATOR 


N = 10 

. 10000  00000  00000 
-.29128  03421  36361 
. 47S73  22804  11881 
-.  36817  31277  22802 
. 1 40?  1 69854  52941 
-.  21678  54057  76489 


00000  +01  , 10000 

44413-01  .12538 

21947  -03  . 762.4? 

06720  -05  .28627 

91735  -07  .67825 

52502  -10  .83325 


00000  00000  00000  +01 
63245  30305  22254  -01 
37336  22270  54117  -04 
47330  74607  14153  -06 
76089  43662  46538  -OS 
19383  71648  11691  -1? 
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TABLE  B. VI 


ERROR  OF  PAPE"  APPROXIMATION  TO  4z~gH(z) 


Let  4z“2Hn(z)  be  the  order  main  diagonal  (see  B.V)  Pade' ap- 
proximant  to  4z_2H(z)  and  define  en(z)  3 |4z"2 llH(z  )-Hn(z ) | . The  tables 
give  en(z)  for  n = 2(2)10  and  z = l(l)lO  . The  expression  in  parentheses 
to  the  right  of  each  number  indicates  the  power  of  10  by  which  the  number  is 
multiplied. 


z\ 

2 

4 

6 

8 

10 

1 

0.584  (-5) 

0.336  (-9) 

0.433  (-14) 

1.000  (-20) 

1.000 

(-20) 

2 

0.443  (-3) 

0.308  (-6) 

0.646  (-10) 

0.535  (-14) 

0.210 

(-18) 

3 

0.414  (-2) 

0.148  (-4) 

0.161  (-  7) 

0.691  (-11) 

0.138 

(-14) 

4 

0.179  (-1) 

0.206  (-3) 

0.728  (-  6) 

0.102  (-  8) 

0.657 

(-12) 

5 

0.599  (-1) 

0.140  (-2) 

0.126  (-  4) 

0.443  (-  7) 

0.721 

(-10) 

6 

0.104 

0.597  (-2) 

0.116  (-  3) 

0.881  (-  6) 

0.308 

(-  8) 

7 

0.177 

0.184  (-1) 

0.680  (-  3) 

0.100  (-  4) 

0.675 

(-  7) 

8 

0.341 

0.441  (-1) 

0.286  (-  2) 

0.751  (-  4) 

0.900 

(-  6) 

9 

0.339 

0.873  (-1) 

0.920  (-  2) 

0.404  (-  3) 

0.812 

(-  5) 

10 

0.410 

0.149 

0.240  (-  1) 

0.167  (-  2) 

0.535 

(-  4) 
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COEFFICIENTS  OF  RATIONAL  APFR QXIMftTI QNS  TO  E(z) 


Here  we  present  coefficients  of  the  numerator  and  denominator  poly- 
nomials of  the  rational  approximation  to  E(z)  defined  in  (5.11).  V/e  give 
coefficients  for  n = 0(l)l0  . The  sequence  of  numbers  given  is  for  the  lowest 
power  to  the  highest  power,  respectively.  The  expression  to  the  right  of  each 
number  is  the  power  of  10  by  which  it  is  multiplied. 

e.g.,  Eg(z)  = +..36. 

z2+12z+36 


NUMERATOR 


N - 00 

. 10000  00000  00000  00000  01 


N - 01 

. 40000  00000  00000  00000  01 

.00000  00000  00000  00000  00 


N = 02 

.36000  00000  00000  00000  02 
.30000  00000  00000  00000  01 
.00000  00000  00000  00000  00 


N a 03 

.48000  00000  00000  00000  03 
.60000  00000  00000  00000  02 
.56666  66666  66666  66667  01 
.00000  00000  00000  00000  00 


DENOMINATOR 


.10000  00000  00000  00000  01 


.4 0000  00000  00000  00000  01 

.10000  00000  00000  00000  01 


.36000  00000  00000  00000  02 

.12000  00000  00000  00000  02 

. 10000  00000  00000  00000  01 


.48000  00000  00000  00000  03 
. 18000  00000  00000  00000  03 
.24000  00000  00000  00000  02 
. 10000  00000  00000  00000  01 


TABIE  B.VII  (Continued) 


NUMERATOR 


N = 04 

.84000  00000  00000  00000  +04 
. 12600  00000  00000  00000  +04 
. 16666  66666  66666  66667  +03 
.41666  66666  66666  66667  +01 
.00000  00000  00000  00000  -99 


- 05 

.18144  00000  00000 
.30240  00000  00000 
.46200  00000  00000 
.21000  00000  00000 
.65666  66666  66666 
.00000  00000  00000 


00000  +06 
00000  +05 
00000  +04 
00000  +03 
66667  +01 
00000  -99 


DENOMINATOR 


.84000  00000  00000  00000  +04 
.33600  00000  00000  00000  +04 
.54000  00000  00000  00000  +03 
.40000  00000  00000  00000  +02 
. 10000  00000  00000  00000  +01 


. 18144 
. 75600 
. 1 3440 
. 1 2600 
. 60000 
. 10000 


00000 

00000 

00000 

00000 

00000 

00000 


00000 

00000 

00000 

00000 

00000 

00000 


00000  +06 
00000  +05 
00000  +05 
00000  +04 
00000  +02 
00000  +01 


N - 06 

. 46569 

60000 

00000 

00000 

+07 

. 46569 

60000 

00000 

00000 

+07 

.83160 

00000 

00000 

00000 

+06 

. 19958 

40000 

00000 

00000 

+07 

. 13776 

00000 

00000 

00000 

+06 

. 37800 

00000 

00000 

00000 

+06 

. 81C00 

00000 

00000 

00000 

+04 

.40320 

00000 

00000 

00000 

+05 

.41 160 

00000 

00000 

00000 

+03 

. 25200 

00000 

00000 

00000 

+04 

. 49000 

00000 

00000 

00000 

+01 

. 84000 

00000 

00000 

00000 

+02 

. 00000 

00000 

00000 

00000 

-99 

. 10000 

00000 

00000 

00000 

+01 

N - 07 

. 13837 

82400 

00000 

00000 

+09 

. 13837 

82400 

00000 

00000 

+09 

. 25945 

52000 

00000 

00000 

+08 

.60540 

48000 

00000 

00000 

+08 

.45276 

00000 

00000 

00000 

+07 

.11975 

04000 

00000 

00000 

+08 

.31416 

00000 

00000 

00000 

+06 

. 13860 

00000 

00000 

00000 

+07 

. 15580 

40000 

00000 

00000 

+05 

. 10080 

00000 

00000 

00000 

+06 

. 46480 

00000 

00000 

00000 

+03 

.45360 

00000 

00000 

00000 

+04 

.71857 

1 4285 

71428 

57143 

+01 

.11200 

00000 

00000 

00000 

+03 

. 00000 

00000 

00000 

00000 

-99 

. 10000 

00000 

00000 

00000 

+01 

60 


n 


TABIE  B.\AII  (Concluded) 


NUMERATOR 


N = 08 

. 46702 
.90810 
. 16432 
. 12612 
.85073 
. 29106 
. 78274 
. 54357 
. 00000 


65600 
72000 
41600 
60000 
60000 
00000 
28571 
1 4285 
00000 


00000 

00000 

00000 

00000 

00000 

00000 

42857 

71428 

00000 


00000  +10 
00000  +05 
00000  +09 
00000  +08 
00000  +06 
00000  +05 
14286  +03 
57143  +01 
00000  -55 


. 17643 
. 35286 
.65585 
. 54054 
.41441 
. 16336 

-Ml 

. 76575 
. 00000 


22560 
45 1 20 
52000 
00000 
40000 
32000 
28571 
85714 
36507 
00000 


00000 

00000 

00000 

00000 

00000 

00000 

42857 

28571 

53650 

00000 


00000  +12 
00000  +11 
00000  +10 
00000  +09 
00000  +08 
00000  +07 
14286  +05 
42857  +03 
75365  +01 
00000  -95 


. 73748 
. 15084 
.28621 
.24872 
. 2021  1 
.50210 
. 36810 
. 76668 
. 12887 
. 58575 
.00000 


68300 

95788 

23264 

04720 

87168 

12000 

65142 

42857 

46031 

36507 

00000 


80000 

80000 

00000 

00000 

00000 

00000 

85714 

14285 

74603 

53650 

00000 


00000  +13 
00000  +13 
00000  +12 
00000  +11 
00000  +10 
00000  +08 
28571  +07 
71429  +05 
17460  +04 
75365  +01 
00000  -99 


X 


DENOMINATOR 


.46702 
.20756 
.42378 
.51891 
.41580 
.22176 
. 75600 
. 14400 
. 10000 


65600 

73600 

33600 

84000 

00000 

00000 

00000 

00000 

00000 


00000 

00000 

00000 

00000 

00000 

00000 

00000 

00000 

00000 


00000  +10 
00000  +10 
00000  +09 
00000  +08 
00000  +07 
00000  +06 
00000  +04 
00000  +03 
00000  +01 


17643 

75394 

16605 

21189 

18162 

10810 

44352 

11880 

18000 

10000 


22560 

51520 

38880 

16800 

14400 

80000 

00000 

00000 

00000 

00000 


00000 

00000 

00000 

00000 

00000 

00000 

00000 

00000 

00000 

00000 


00000  +12 
00000  +11 
00000  +11 
00000  +10 
00000  +09 
00000  +08 
00000  +06 
00000  +05 
00000  +03 
00000  +01 


. 73748 
. 33522 
.71455 
.94057 
.84756 
. 54486 
.25225 
.82368 
. 17820 
. 22000 
. 10000 


68300 
1 2864 
06366 
20320 
67200 
43200 
20000 
00000 
00000 
00000 
00000 


80000 

00000 

00000 

00000 

00000 

00000 

00000 

00000 

00000 

00000 

00000 


00000  +13 
00000  +13 
00000  +12 
00000  +11 
00000  +10 
00000  +09 
00000  +08 
00000  +06 
00000  +05 
00000  +03 
00000  +01 


TABUS  B. VIII 


r 


ERRCK  CF  RATIONAL  APPROXIMATION  TO  E(z) 

Let  En(z ) be  the  order  rational  approximation  to  E(z)  as 
defined  in  (5.1l)  and  Table  B.VII.  Set  en(z)  = |E(z)-En(z)  | . The  tables 
give  en(z)  for  n = 3(l)lO  , r = l(l)lO  and  0 = 0,n/2^n  where  z = rei0  . 
The  expression  in  parentheses  to  the  right  of  each  number  indicates  the  power 
of  10  by  which  the  number  is  multiplied. 


r\ 

0 

n/2 

n 

n (Relative  Error) 

n = 

1 

0.593  (-5) 

0.177  (-4) 

0.487  (-4) 

0.369  (-4) 

2 

0.272  (-4) 

0.291  (-3) 

0.219  (-2) 

0.119  (-2) 

3 

0.134  (-4) 

0.148  (-2) 

0.317  (-1) 

0.115  (-1) 

0.870  (-4) 

0.444  (-2) 

0.284 

0.643  (-1) 

5 

0.286  (-3) 

0.975  (-2) 

0.175  ( 1) 

0.230 

6 

0.566  (-3) 

0.172  (-1) 

0.719  ( 1) 

0.517 

7 

0.895  (-3) 

0.259  (-1) 

0.206  ( 2) 

0.762 

0 

0.125  (-2) 

0.345  (-1) 

0.490  ( 2) 

0.896 

9 

. 0.160  (-2) 

0.413  (-1) 

0.110  ( 3) 

0.957 

10 

0.193  (-2) 

A 

0.454  (-1) 

0.243  ( 3) 

0.976 

• n = 

i 

Ht 

0.108  (-6) 

0.287  (-6) 

0.754  (-6) 

0.571  (-6) 

2 

0.407  (-5) 

0.543  (-5) 

0.703  (-4) 

0.382  (-4) 

3 

0.226  (-4) 

0.155  (-4) 

0.180  (-2) 

0.655  (-3) 

4 

0.716  (-4) 

0.198  (-3) 

0.276  (-1) 

0.624  (-2) 

5 

0.166  (-3) 

0.102  (-2) 

0.326 

0.429  (-1) 

6 

0.315  (-3) 

0.317  (-2) 

0.374  ( 1) 

0.269 

7 

0.524  (-3) 

0.724  (-2) 

0.274  ( 3) 

0.101  ( 2) 

8 

0.787  (-3) 

0.133  (-1) 

0.773  ( 2) 

0.141  ( 1) 

9 

0.110  (-2) 

0.209  (-1) 

0.126  ( 3) 

0.110  ( 1) 

10 

0.145  (-2) 

0.288  (-1) 

0.256  ( 3) 

0.103  ( 1) 

TAB IE  B. VIII  (Continued) 


l 


► 


N? 

£_ 

0 

nil 

TT 

n (Relative  Error) 

n = 

_5  0 

1 

0.376  (-8) 

0.632  (-8) 

0.125  (-7) 

0.947 

(-8) 

2 

0.140  (-6) 

0.379  (-6) 

0.188  (-5) 

0.102 

(-5) 

3 

0.918  (-6) 

0.449  (-5) 

0.674  (-4) 

0.245 

(-4) 

4 

0.290  (-5) 

0.3C6  (-4) 

0.148  (-2) 

0.335 

(-3) 

5 

0.598  (-5) 

0.145  (-3) 

0.242  (-1) 

0.318 

(-2) 

6 

0.916  (-5) 

0.5 01  (-3) 

0.314 

0.226 

(-D 

7 

0.107  (-4) 

0.136  (-2) 

0.318  ( 1) 

0.118 

8 

0.860  (-5) 

0.298  (-2) 

0.216  ( 2) 

0.395 

9 

0.117  (-5) 

0.553  (-2) 

0.848  ( 2) 

0.737 

10 

0.129  (-4) 

0.895  (-2) 

0.228  ( 3) 

0.916 

n = 6 


1 

0.760  (-10) 

0.125  (-9) 

0.225  (-9) 

0.170  (-9) 

2 

0.6C2  (-8) 

0.154  (-7) 

0.559  (-7) 

0.304  (-7) 

3 

0.665  (-7) 

0.255  (-6) 

0.233  (-5) 

0.847  (-6) 

4 

0.337  (-6) 

0.184  (-5) 

0.603  (-4) 

0.136  (-4) 

5 

0.113  (-5) 

0.883  (-5) 

0.129  (-2) 

0.170  (-3) 

6 

0.295  (-5) 

0.290  (-4) 

0.229  (-1) 

0.165  (-2) 

7 

0.649  (-5) 

0.936  (-4) 

0.343 

0.127  (-1) 

8 

0.126  (-4) 

0.288  (-3) 

0.465  ( 1) 

0.850 

9 

0.221  (-4) 

0.778  (-3) 

0.825  ( 2) 

0.717 

10 

0.360  (-4) 

0.178  (-2) 

0.508  ( 3) 

0.204 

n = 7 


1 

0.144 

(-11) 

0.236  ( 

-11) 

0.407 

(-11) 

0.308 

(-11) 

2 

0.229 

(- 

9) 

0.579  ( 

- 9) 

0.185 

(- 

8) 

0.101 

(- 

e) 

3 

0.370 

(- 

8) 

0.137  ( 

- 7) 

0.936 

(- 

7) 

0.340 

(- 

7) 

4 

0.236 

(- 

7) 

0.113  ( 

- 6) 

0.240 

(- 

5) 

0.543 

(- 

6) 

5 

0. 905 

(- 

7) 

0.526  ( 

- 6) 

0.551 

(- 

4) 

0.725 

(- 

5) 

6 

0.251 

(- 

6) 

0.142  ( 

- 5) 

0.120 

(- 

2) 

0.8S3 

(- 

4) 

7 

0.553 

(- 

6) 

0.627  ( 

- 5) 

0.229 

(- 

1) 

0.848 

(- 

3) 

8 

0.1C5 

(- 

5) 

0.341  ( 

- 4) 

0.367 

0.671 

(- 

2) 

9 

0.172 

(- 

5) 

0.130  ( 

- 3) 

0.493 

( 

1) 

0.429 

(- 

1) 

10 

0.255 

(- 

5) 

0.380  ( 

- 3) 

0.506 

( 

2) 

0.203 
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r\ 

TABIE 

B.VIII  (Concluded) 

K 

E_ 

n = 

0 

_8 

v/2 

TT 

tt  (Relative  Error 

1 

0.259  (-13) 

0.420 

(-13) 

0.712 

(-13) 

0.539 

(-13) 

2 

0.825  (-11) 

0.208 

(-10) 

0.626 

(-10) 

0.340 

(-10) 

3 

0.201  (-  9) 

0.759 

(-  9) 

0.431 

(-  8) 

0.157 

(-  8) 

4 

0.173  (-  8) 

0.939 

(-  8) 

0.116 

(-  6) 

0.262 

(-  7) 

5 

0.853  (-  8) 

0.648 

(-  7) 

0.237 

(-  5) 

0.312 

(-  6) 

6 

0.297  (-  7) 

0.327 

(-  6) 

0.528 

(-  4) 

0.380 

(-  5) 

7 

0.820  (-  7) 

0.142 

(-  5) 

0.120 

(-  2) 

0.444 

(-  4) 

8 

0.192  (-6) 

0.555 

(-  5) 

0.242 

(-  1) 

0.442 

(-  2) 

9 

0.399  (-  6) 

0.189 

(-  4) 

0.418 

0.363 

(-  2) 

10 

0.754  (-  6) 

Q 

0.558 

(-  4) 

0.638 

( 1) 

0.256 

(-  1) 

n = 

i 

0.439  (-15) 

0.709 

(-15) 

0.119 

(-14) 

0.902 

(-15) 

2 

0.280  (-12) 

0.708 

(-12) 

0.2C6 

(-11) 

0.112 

(-11) 

3 

0.102  (-10) 

0.391 

(-10) 

0.205 

(-  9) 

0.745 

(-10) 

4 

0.117  (-  9) 

0.655 

(-  9) 

0.657 

(-  8) 

0.149 

(-  8) 

5 

0.710  (-  9) 

0.571 

(-  8) 

0.127 

(-  6) 

0.167 

(-  7) 

6 

0.290  (-8) 

0.330 

(-  7) 

0.235 

(-  5) 

0.169 

(-  6) 

7 

0.904  (-  8) 

0.142 

(-  6) 

0.535 

(-  4) 

0.198 

(-  5) 

8 

0.231  (-  7) 

0.474 

(-  6) 

0.126 

(-  2) 

0.230 

(-  4) 

9 

0.509  (-7) 

0.124 

(-  5) 

0.266 

(-  1) 

0.231 

(-  3) 

10 

n = 

0.995  (-  7) 

10 

0.283 

(-  5) 

0.484 

0.194 

(-  2) 

1 

0.658  (-17) 

0.114 

(-16) 

0.188 

(-16) 

0.142 

(-16) 

2 

0.898  (-14) 

0.227 

(-13) 

0.646 

(-13) 

0.351 

(-13) 

3 

0.492  (-12) 

0.189 

(-H) 

0.952 

(-11) 

0.346 

(-11) 

4 

0.747  (-11) 

0.425 

(-10) 

0.392 

(-  9) 

0.887 

(-10) 

5 

0.566  (-10) 

0.463 

(-  9) 

0.831 

(-  8) 

0.109 

(-  8) 

6 

0.267  (-  9) 

0.315 

(-  8) 

0.133 

(-  6) 

0.956 

(-  8) 

7 

0.101  (-8) 

0.151 

(-  7) 

0.241 

(-  5) 

0.893 

(-  7) 

8 

0.296  (-8) 

0.546 

(-  7) 

0.569 

(-  4) 

0.104 

(-  5) 

9 

0.742  (-  8) 

0.154 

(-  6) 

0.139 

(-  2) 

0.121 

(-  4) 

10 

0.164  (-  7) 

0.426 

(-  6) 

0.256 

(-  1) 

0.103 

(-  3) 
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APPENDIX  C 


FORTRAN  PROGRAM FCR  COMPUTING  RATIONAL  APPROXIMATIONS 


Here  we  describe  a FORTRAN  program  which  computes  the  rational  ap- 
proximations (2.1)  and  (2.18)  for  the  incomplete  gamma  function  and  its  spe- 
cial cases.  The  selection  of  input  data  determines  the  function  to  be  approxi- 
mated. We  also  include  a description  of  input  data,  operating  procedures, 
output,  and  a listing  of  the  FORTRAN  program. 

The  input  data  are  read  in  the  order  r,  0,  A,  B,  NCODE,  IOPT,  LOFT, 


1“  1 
L Error  J 


r n i 

The  value  IOPT  determines  which  value  in  the  braces  < or  s 

l Error  J 

should  be  read.  If  the  ntb  approximant  is  desired,  IOFT  * 7 is  entered  and 
the  corresponding  value  in  braces,  n , is  entered.  IOPT  / 7 instructs  the 
computer  to  iterate  until  an  integer  m is  found  6uch  that 

|v  zuez  ||Vm(z,v)  - Vm  _l(z,v)|  < Error  , and  therefore  Error  should  accompany 
this  choice  of  IOPT.  Selection  of  LOFT  offers  the  choice  of  computing  either 
vn(z>u)  or  sn(z>v)  (see  (2.1)  and  (2.18)).  If  LOP T = 9 , the  rational 
approximate  is  Vn(z,v)  and  if  LOFT  £ 9 , the  approximate  is  Sn(z,v)  . 

For  z = re1®  (£  in  degrees)  and  v = A+iB  , the  following  table 
indicates  the  values  of  u and  NCODE  to  select  to  compute  the  approximations 
to  the  designated  functions.  A listing  of  the  FORTRAN  program  concludes  this 
Appendix. 
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il  '5 

31  * 


•s 

C4 

>— * 

OJ 

O 

♦“3 

• 

1 

rH 

HCY 

►“3 

^ 1 

CVJ 

CVJ 

p 

hIc\j 

' — 1 

- — 1 

II 

II 

II 

-p 

■p 

+> 

t* 

-p 

•p 

■p 

a 

(0 

a 

•H 

o 

•H 

.CO 

o 

U) 

i|(\l 

^KM 

■ 

r^M 

■P 

■ 

•p 

-P 

ml  m 

<1  < 


{ vn^z2e'irr^2^)+Rn(z2e'in/2,i)} 


c 


203 

101 


102 


103 

90 


80 

81 


70 


1 


100 


105 


110 


RATIONAL  APPROX. ,A$CEN0 I NG  SERIES 
DIMENSION  VR(40),VI (40) 

READ  400,R,TH,A,B,NC0DE, I OPT.KOPT, LOPT, 
PRINT  403,R.TH, A,B 
Pl=3. 1415926 
THA=TH*PI/180. 

HPI=.5*PI 
RTPI=1. 7724538 
NN=  1 

X=R*C0SF(THA) 

Y=R*S1NF(THA) 

RO=LOGF(R) 

IF  (IOPT-7)  101,102,101 

READ  401, ERROR 

ERRCR=ERROR*ERROR 

GO  TO  103 

READ  402, NT 

MT=NT+1 

GO  TO  (70,80,50, 100,1 10), NCODE 
CM=.  5 
GO  TO  81 
CM= 1 . 0 

RZ=(Y*Y-X*X)*CM 
RZI=-2.*X*Y*CM 
THMTHA+R7. 1 
C=COSF(TH) 

D=S I NF(TH) 

DEN=R*EXPF(RZ) 

CR=C*DEN 

CI  = D*DEN 

X=RZ 

Y=R7  I 

GO  TO  300 

C=»A*RO-B*THA+X 

DN= 1 . / ( A*A+B*B ) 

TH=B*RO+A*THA+Y 
DN=EXPF (C)*DN 

CR=DN*(A*COSF(TH)+B*S I NF(TH) ) 
CI»DN*(A*S I NF(TH)-B*COSF(TH) ) 

GO  TO  300 

THA-.  5*THA 

OX 

X=-Y 

Y=C 

DN«SQRTF(R)*EXPF(X) 

CR“ DN*COSF ( THA+Y ) 

CI*DN*SINF(THA+Y) 

GO  TO  300 
OX*X-Y*Y 
D=-2.  *X*Y 
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I 


Y=C 

X-D 

106  DN-R*EXPF( X) 

CR=DN*C0SF(THA+Y ) 
CI=DH*SINF(THA+Y) 

GO  TO  300 

60  IF  (NN-2)  61,62,6?. 

61  NN-2 
X— X 
Y— Y 
CC-C 
DD-D 

IF(NC0DE-*:)106,1O5.106 

300  IF  (LOPT-i)  202,205,202 

202  A-A-1 . 

205  AI»(A+1.  )*(A+1.  )+B*B 

AI-1.  / A I 

I F(L0PT-5 ) 1 1 , 10, 1 1 

10  PI-1. 

S2-A+2.  +X 
SI-1. 

P2-(A+2.  )-( ( A+1 . )*X+B*Y)*AI 
Q1“0. 

Q2-B+(B*X-Y*(A+1.  ))*AI 

T 1-0. 

T2-D+Y 
GO  TO  301 

11  PI-0.0 
P2-A+2. 

OJ-O.O 

Q2-D 

51- 1.0 

52- X+A+2, 

T 1-0.  0 
T2-Y+B 

301  J=1 

VR(  1 )=(P2*S2+Q2*T2 )/ ( S2*S2+T2*T2 ) 
VI (1 )=(Q2*S2-P2*T2)/(S2*S2+T2*T2) 

302  K-J+1 
F-J 
G-K 

H-2.  *F+A 

AL=M*(G+A)-B*B 

BL=B*(H+G+A) 

C - ( H+ 1 . )*(H*(H+2. )-B*B+A*X-B*Y) 
C=C-B*(B*X+A*Y+B*(2.*H+2.  )) 
D-B*(H*(H+2.  )-B*B+A*X-B*Y) 
D-D+(H+1 . )*(0*X+A*Y+B*(2.  *H+2.  )) 
E-F*( (X*X-Y*Y)*(H+2,  )-2.*B*X*Y) 
FE- F* ( 2.  *X*Y*{ H+2.  )+B* ( X*X-Y*Y ) ) 
DI4—  1 . / ( AL*AL+BL*BL ) 
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303 

305 

304 

306 

307 

308 

305 

310 

30 

40 


50 


ACPBD=AL*C+BL*D 

BCMAD=BL*C-AL*D 

AEPBF=AL*E+BL*FE 

BEMAF=BL*E-AL*FE 

P3=Dil*  ( P2*ACPBD+Q2*BCMAD+P  1 *AE  PBF+Q 1 *G  EMAF ) 

Q3=DN*  (-P2*BCMAD+Q2*AC  PBD-P 1 *BEHAF+Q1 *AEPBF ) 

S3=Di'l*(S2*ACPBD+T2*BCMAD+S1*AEPBF+T1*BEMAF) 

T3«DN*(-S2*BCMAD+T2*ACPBD-S 1 *BEMAF+T 1 *AEPBF) 

Q1=Q2 

Q2»03 

P1  = P2 

P2=P3 

51- S2 

52- S3 
T1-T2 
T2-T3 
L-K+1 

IF  (I0PT-7)  304,303,304 
IF  (NT-L)  304,304,305 
J«J+1 
GO  TO  302 

DEN-1. /(S2*S2+T2*T2) 

VR(K)-(P2*S2+Q2*T2)*DEN 
VI (K)-(Q2*S2-P2*T2 )*DEN 
IF  (IOPT-7)  300,307,308 
IF  (NT-K)  308,308,305 
RER-VR(K)-VR(K-1 ) 

RE  I -VI (K)-VI (K-1 ) 

IF  (IOPT-7)  309,310,305 
C-CR*RER-CI*REI 
D»CI*RER+CR*REI 
ERT-C*C+D*D 

IF  (ERT-ERROR)  310,310,305 
C»CR*VR(K)-CI*VI (K) 

D-CI*VR(K)+CR*VI (K) 

GO  TO  (30,40,50,60,60),NC0DE 
PRINT  500, 

PRINT  501,K,C,D 
GO  TO  68 
PRINT  502, 

PRINT  501,K,C,D 
C-. 88622692-C 
0—0 

PRINT  501 ,K, C, 0 
GO  TO  68 
PRINT  502, 

PRINT  501 ,K, C, D 
C-1.  2533141-C 
D— I) 

PRINT  501 ,K,C, D 
GO  TO  68 


70 


62 

64 

65 


67 

66 


60 

65 

400 

401 

402 

403 

500 

501 

502 

505 

506 


CM=  1 . 

IF  (NCODE-4)  64,65,64 
CM=.  5 

SIZR*CM*(DD-D) 

SIZ|»CM*(C-CC) 

CIZR=CI1*(C+CC) 

CIZI=CM*(D+DD) 

PRINT  505, 

PRINT  506,K,SIZR,5IZ|,CIZR, CIZI 

DN»1. 2533141 

IF  (NC0DE-5)  66,67,66 

DN-.  5*DN 

CIZR=DN-CIZR 

S I ZR=*DN-S  I ZR 

CIZI =-C I Z I 

SIZI=-SIZI 

PR I NT  506, K, S I ZR, S I Z I , C I ZR, CIZI 
IF  (SENSE  SWITCH  1)  65,1 
PRINT  501,K,VR(K),VI (K) 

GO  TO  1 

FORMAT (4F 12. 0,4| 3) 

FORMAT ( E14. 7 ) 

FORMAT ( I 3) 

FORMAT (4( E 14. 7 , 2 X ) / ) 

FORMAT (15H  INC.  GAMMA  FM ./) 
FORMAT ( I 3, 1 X,2(E14. 7,2X)/) 
FORMAT ( 15H  ERROR  FUNCTION/) 
FORMAT (1011  FRESNEL  INTEGRALS/) 
FORMAT ( I 3, 2X,4(E14. 7.2X)/) 

END 


APPENDIX  D 


FORTRAN  PROGRAM  FOR  COMPUTATION  OF  THE  EXPONENTIAL 


AND  CIRCULAR  FUNCTIONS 


The  following  FORTRAN  program  is  based  on  the  results  of  Section  IV. 

Input : (Cards)  l£t  n be  the  order  of  the  rational  approximations 
to  the  exponential  and  circular  functions,  and  suppose  that  these  approxima- 
tions are  desired  for  n * ni(l  )nk  and  x * x^(^x)xm  . The  data  are  entered 
in  the  order  nx  > nk  , x1  , , £x  . 

Switch  settings : None . 

Qitput : (Printed)  The  output  is  of  the  following  form  for 

xi>X2,  •••  >x jj. 


en11  > sinn1(Jtl)  > cofin1^xl)  > tann1(xl> 
en^  > 6lnn2(xl)  » cosn2^xl^  > tann2(xl) 


(x  ) 

enk1  > 6lnnk(xl)  > cosnk(xl)  > ta%(xl) 


c 

1 


2 


3 


5 

6 

100 

200 

300 


RATIONAL  APPROX.  TO  EXP(-Z ) , COS(Z ), S I N(Z ), TAN(Z ) 

READ  100.NI ,NF,Z| ,ZF,ZD 

PRINT  200 

7=7.  I 

ZP»Z*Z 

ZM— ZP 

N=2 


FM=2.0 

CM-2.*(2.*FM+1.  ) 

A1P-2.  0 

A2P«1 2.  O+ZP 

B 1 P=  1 , 0 

B2P=6„  0 

A1M=2,0 

A2H-12,  0+7.H 

B1M»1k0 

B2M=6.  0 

N®N+1 

A3P-CM*A2P+ZP*A1P 

A3M»CM*A2M+ZM*A1M 

B3P=CM*B2P+ZP*B1P 

B 3M®  CM*B  2M+Z  M*  B 1 M 

A1P-A2P 

A2P®A3P 

B1P-B2P 

B2P»B3P 

A1H-A7.M 


A2M»A3M 

B1M-B2M 

B2M-B3M 


FM=FM+1t  0 
CM=2. *(2„  *FM+1 . 0) 

IF  (Nl-N)  4,4,3 

EZ=(A3P-Z*B3P)/(A3P+7.*B3P) 

D=A3M*A3M+Z  P*B3H*B3M 

CZ= ( A3M*A3M+ZM*B3M*B3M ) / D 

SZ=2.*Z*A3M*B3M/D 

D=S7./CZ 

PRINT  300,N,7,EZ,SZ,CZ,D 
IF  (NF-N)  1,5,3  - 
IF  (7.F-Z)  1,1,6 
Z-7+ZD 


GO  TO  2 

FORMAT (2(I3),3(E15. 7)) 

F0RMAT(//2H  N, 7X1  HZ,  1 3X7HEXP(-Z  ),9X6HS  I N(7. ),  1 0X6HC0S(7 ) 
1, 10X6HTAN(Z)/) 

FORMAT ( I 2,5(2XE14. 7) ) 

END 


APPENDIX  E 


FORTRAN  PROGRAM  FOR  COMPUTING  RATIONAL  APPROXIMATIONS  TO  E(z 


Here  we  give  a brief  description  of  a FORTRAN  program  which  computes 
the  rational  approximations  to  E(z)  defined  by  (5.11).  A listing  of  the 
program  follows  the  description  of  input  and  output  data. 

Let  z = re1-®  and  suppose  we  wish  to  compute  the  ntb  order  rational 
approximation  to  E(z)  as  defined  by  (5.11)  for  rj(Ar)rF  and  0-j-(A0)0p 
where  I and  F subscripts  denote  initial  and  final  values,  respectively. 

Input : (Cards ) The  values  r-j-  , rF  , Ar  , 0j  , 0F  and  AO  are 
entered  in  this  order  (0  in  degrees)  , three  numbers  per  card. 

Output:  (Typed)  For  each  pair  of  values  of  r and  9 , the  n^*1 


order  rational  approximation  is  printed,  real  part  and  then  the  imaginary 
part. 


c rational  approximations  to  e(z) 

DIMENSION  F1(30)tF2(30),F3(30),F4(30),P1(30), 
P2(30),P3(30),P4(30) 

1 READ  100.RI ,RF,RD,THItTHF,THD 
CF-3. 141592653589793238462643/100. 

DO  2 J«1,30 

F1(J)-0.0 
F2( J)“0. 0 
F3(J)-0.0 
F4(J)»0„  0 
P1(J)-0.0 
P2(J)-0.0 
P3(J)=0.0 

2 P4(J)»0.0 
F3(O=1.0 
F3(2)*24.0 
F3(3)-180.0 
F3(4)-480.0 
F2(1)=  1.0 
F2(2)-12.0 
F2(3)=36.0 
FI ( 1 )®1„ 0 
FI ( 2 )=4. 0 
P3  ( 1 )=0.  0 
P3(?)=17. /3. 

P3(3)»60.0 
P3(4)=480. 0 
P2( 1 )=0„ 0 
P2(2)*3.0 
P2(3)»36.0 
P1(1)=0.0 
P1(2)»4.0 
N=4 

3 M*=N+I 
H=N 

TWH=2.*H 

A1=(H-2.  ) * ( T V/H—  1 «,  )/(H*(TWH-3.  )) 

A2-1.0 
A3— A1 

B1  = 2.*(TWH-1.  )*(H+1. )/H 

02— 2.*(TWH-1.  )*(H-3.  )/H 

F4( 1 )-A1*F3( 1 )+A2*F2( 1 )+A3*F1 ( 1 ) 

P4 ( 1 )«A 1 *P3 ( 1 )+A2*P2 ( 1 ) +A3*P 1 ( 1 ) 

DO  4 J»2,M 

F4(J)»AI*F3(J)+B1*F3(J-1 )+A2*F2(J)+B2*F2(J-1 )+A3*F1 (J) 

4 P4(J)»A1*P3(J)+B1*P3(J-1 )+A2*P2( J)+B2*P2( J-1 )+A3*P1 (J) 
DO  5 J"1 tM 

FI (J)»F2( J) 

F2(J)-F3(J) 

F3(J)-F4(J) 

P1(J)-P?.(J) 


75 


1 


P2(J)=P3( J) 

5 P3(J)-P4(J) 

6 RO-RI 

7 R=1./RO 
TH-THI 

8 PRINT  101,R0,TH 
T-TH*CF 
RE»R*COSF(T) 

RIM— R*SINF(T) 
SRN«P4(M)*RE+P4(M-1 ) 
SIN«P4(M)*RIM 
SRD«F4(M)*RE+F4(M-1 ) 
SID-F4(M)*RIM 

DO  S’  J=3,M 
L=M-J+1 

S»RE*SRN-R I M*S I N+P4 ( L ) 

S I N=RE*S I N+R I M*SRN 
SRN=S 

S«RE*SRD-RIM*SID+F4(|_) 

S I D=RE*S I D+R I M*SRD 
? SRD=S 

DEN=SRD*SRD+S I D*S I D 
QR= ( SRN*SRD+S I N*S I D ) /DEN 
Q I = ( S I N*SRD-SRN*SI D)/DEN 
PRINT  1 04,QR,Q I 
TH=TH+THD 
IF  (THF-TH) 10,8,8 
10  RO-RO+RD 

IF  (RF-RO) 1 2,7,7 
17.  CONTINUE 
13  N=N+1 

GO  TO  3 

100  FORMAT  (3E15.7) 

101  FORMAT ( El 4,  7, 2 X, El 4. 7 ) 
104  FORMAT(2(E3?..  24)) 

END 


r 


; 


DISTRIBUTION  LIST 
Contract  No.  Nonr-2638(00)(x ) 


No.  of  Copies  Recipient 

10  Chief,  Bureau  of  Ships 

Washington  25,  D.  C. 

3 Tech.  Info  Br  (Code  335) 

1 Tech.  Asst,  to  Chief  (Code  106) 

3 Computer  Systems  and  Applications  Br 
(Code  732) 

1 Asst.  Chief  for  Field  Activities 
(Code  700) 

1 Asst.  Chief  for  Nuclear  Propulsion 
(Code  1500) 

1 Lab  Mgt.  (Code  320) 

1 Chief,  Bureau  of  Weapons 

18th  and  Constitution  Avenue 
Washington  25,  D.  C. 

1 Chief,  Bureau  of  Supplies  and  Accounts 

18th  and  Constitution  Avenue 
Washington  25,  D.  C. 

1 Chief,  Bureau  of  Census 

Washington  25,  D.  C. 


1 


Chief,  Office  of  Naval  Research 
Washington  25,  D.  C. 


1 Naval  Shipyard 

Boston,  Massachusetts 

1 Naval  Shipyard 

Charleston,  South  Carolina 

1 Naval  Shipyard 

Long  Beach,  California 

1 Naval  Shipyard 

Attn:  Shipyard  Tech.  Library 
Bldg.  746,  Code  370a 
Mare  Island,  California 


. 


1 


I 


h 


DISTRIBUTION  LIST  (Continued) 

No.  of  Copies  Recipient 

1 Naval  Shipyard 

New  York,  New  York 

1 Naval  Shipyard 

Norfolk,  Virginia 

1 Naval  Shipyard 

San  Francisco,  California 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Naval  Shipyard 
Philadelphia,  Pennsylvania 

Naval  Shipyard 
Portsmouth,  New  Hampshire 

Naval  Shipyard 

Puget  Sound,  Washington 

Naval  Shipyard 
Pearl  Harbor,  Hawaii 

CO  & DIR,  USNBTL 
Philadelphia,  Pennsylvania 

CO  & DIR,  U.  3.  Naval  Electronics  Lab 
San  Diego,  California 

CO  & DIR,  U.  S.  Naval  Research  and 
Development  Lab 
San  Francisco,  California 

CO  & DIR,  U.  S.  Naval  Training  Device 
Center 

Computer  Branch 

Port  Washington,  New  York 

CO,  USNCML 

St.  Paul,  Minnesota 

CDR,  U.  S.  Naval  Proving  Ground 
Dahlgren,  Virginia 


2 


No.  of  Copies 
3 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


3 


DISTRIBUTION  LI3T  (Continued) 


Recipient 

CDR,  U.  S.  Naval  Ordnance  Test  Station 
China  Lake,  California 
1 Pasadena  Annex 
1 Tech.  Library 
1 Michelson  Lab  (Code  5038) 

CDR,  U.  S.  Naval  Ordnance  Lab 
White  Oak,  Silver  Spring,  Maryland 

Director,  United  States  Naval  Academy 
Annapolis,  Maryland 

DIR,  U.  S.  Naval  Research  Lab 
18th  and  Constitution  Avenue 
Washington  25,  D.  C. 

SUPT,  U.  3.  Naval  Postgraduate  School 
Monterey,  California 

Attn:  Library,  Technical  Reports  Section 

SUPT,  U.  S.  Naval  Academy 
Department  of  Math 
Annapolis,  Maryland 

CG,  Aberdeen  Proving  Ground 
Aberdeen,  Maryland 

DIR,  National  Bureau  of  Standards 
Washington  25,  D.  C. 

G.  C.  Marshall  Flight  Center 
NASA 

Huntsville,  Alabama 

DIR,  Langley  Aero  Lab 
NASA 

Langley  Field,  Virginia 

DIR,  Lewis  Flight  Propulsion  Lab 
NASA 

Cleveland  11,  Ohio 


— 


* ■ < 


I 


DISTRIBUTION  LIST  (Continued) 

Recipient 

Commander 

Aeronautical  Systems  Division 
Attn:  Applied  Mathematics  Laboratory 
Wright-Patterson  Air  Force  Base,  Ohio 

CG,  Frankford  Arsenal 
Head,  Math  Section 
Pitman-Dunn  Laboratories 
Philadelphia  37,  Pennsylvania 

CG,  White  Sands  Proving  Ground 
Flight  Determination  Lab 
Las  Cruces,  New  Mexico 

USAEC,  Technical  Information  Service 
Oak  Ridge,  Tennessee 

U.  3.  Atomic  Energy  Commission 
Washington  25,  D.  C. 

Attn : Tech . Library 

CO,  Diamond  Ordnance  Fuze  Laboratory 
Washington  25,  D.  C. 

Attn : Library 

George  Washington  University 
Logistics  Research 
Washington,  D.  C. 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
Silver  Spring,  Maryland 

University  of  California 
Librarian,  Numerical  Analysis 
Los  Angeles  24,  California 

Carnegie  Institute  of  Technology 
Pittsburgh,  Pennsylvania 


DISTRIBUTION  LIST  (Continued) 


Recipient 


Hudson  Laboratory 
Columbia  University 
Dobbs  Ferry,  New  York 

Harvard  University 
Cambridge,  Massachusetts 
1 Department  of  Math 

Attn:  Prof.  J.  S.  Wal6h 
1 Computation  Laboratory 

Institute  for  Advanced  Study 
Princeton,  New  Jersey 

Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 
Attn:  Computation  Center 

New  York  University 
New  York,  New  York 

1 Inst,  of  Math  Sciences 
1 AEC  Computing  Facility 

DIR,  Research  Center 
Ohio  State  University 
Columbus,  Ohio 

Pennsylvania  State  University 
Department  of  Math 
University  Park,  Pennsylvania 

Princeton  University 
Princeton,  New  Jersey 
Attn : Library 

University  of  California 
Berkeley,  California 

University  of  Illinois 
Urbana,  Illinois 

1 Department  of  Math 
1 Electronic  Digital  Computer  Lab 


DISmiEUTION  LIST  (Continued) 


f * 


No.  of  Copies 
1 


Recipient 

University  of  Kansas 
Department  of  Mathematics 
Lawrence,  Kansas 

University  of  Maryland 
College  Park,  Maryland 
1 Department  of  Math 
1 Inst,  for  Fluid  Dynamics  and 
Applied  Math 

1 Attn:  Professor  Burt  Hubbard 

Willow  Run  Research  Center 
University  of  Michigan 
Ypsilanti,  Michigan 

University  of  Michigan 
Ann  Arbor,  Michigan 

University  of  Missouri 
Department  of  Mathematics 
Columbia,  Missouri 

University  of  Washington 
Department  of  Mathematics 
Seattle,  Washington 

Yale  University 

New  Haven,  Connecticut 

State  College  of  Washington 
Department  of  Mathematics 
Pullman,  Washington 

Stanford  University 
Palo  Alto,  California 
Attn:  Dr.  G.  E.  Forsythe 

Johns  Hopkins  University 
Charles  and  34th  Streets 
Baltimore  18,  Maryland 


DISTRIBUTION  LIST  (Continued) 


h 


No.  of  Copies 
1 


Recipient 


Rutgers  University 

New  Brunswick,  New  Jersey 

Attn:  Professor  E.  P.  Starke 

Brown  University 
Division  of  Engineering 
Providence,  Rhode  Island 
Attn:  Dr.  R.  D.  Kodis 

University  of  Rochester 
Department  of  Mathematics 
Rochester,  New  York 

The  Scripps  Institute  of  Oceanography 
University  of  California 
LaJolla, California 
Attn:  Dr.  Walter  Munk 

Dr.  Werner  C.  Rheinboldt,  Director 
Computing  Center,  112  Hinds  Hall 
Syracuse  University 
Syracuse  10,  New  York 

DIR,  Combustion  Engineering,  Inc. 
Windsor,  Connecticut 

DIR,  Westinghouse  Electric  Corporation 
Bettis  Atonic  Power  Division 
P.  0.  Box  1468 
Pittsburgh  30,  Pennsylvania 

Argonne  National  Laboratory 
P.  0.  Box  299 
Lemont,  Illinois 

Armour  Research  Foundation 
35  West  33rd  Street 
Chicago  16,  Illinois 

Battelle  Memorial  Institute 
505  King  Avenue 
Columbus,  Ohio 


DISTRIBUTION  LIST  (Continued) 


1 


Brookhaven  National  Laboratory 
Upton,  Long  Island,  New  York 


1 


Cornell  Aero  Laboratory,  Inc . 
Buffalo,  New  York 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Douglas  Aircraft  Company,  Inc . 

Santa  Monica  Division 
3000  Ocean  Park  Boulevard 
Santa  Monica,  California 

Engineering  Research  Associates 
St.  Paul  4,  Minnesota 

IBM  Corporation 
590  Madison  Avenue 
New  York  22,  New  York 

Knolls  Atomic  Power  Laboratory 
General  Electric  Company 
Math  Analysis  Unit 
Schenectady,  New  York 

Lincoln  Lab 
B-125 

Lexington,  Massachusetts 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division 
Sunnyvale,  California 

Lockheed  Aircraft  Corporation 
Van  Nuys,  California 

Los  Alamos  Scientific  Lab 
Los  Alamos,  New  Mexico 

Remington  Rand  Division  of  Sperry  Rand 
Electronic  Computer  Department 
315  Fourth  Avenue 
New  York  10,  New  York 


8 


I 


DISTRIBUTION  LIST  (Continued) 


M — 


No.  of  Copies 


Recipient 


1 


1 


1 


1 


1 


1 


2 


1 


1 


Ramo  Wooldridge  Corporation,  a Division 
of  Thompson  Ramo-Wooldridge,  Inc. 
8433  Fallbrook  Avenue 
Canoga  Park,  California 

Rand  Corporation 
Santa  Monica,  California 

Sandia  Corporation 
Albuquerque,  New  Mexico 
Attn : Library 

United  Aircraft  Corporation 

400  Main  Street 

East  Hartford  8,  Connecticut 

Division  of  Vitro  Labs 
Vitro  Corporation  of  America 
Silver  Spring,  Maryland 

Oregon  State  College 

Covallis,  Oregon 

Attn:  Professor  W.  Milne 

Mathematics  Research  Center 
University  of  Wisconsin 
Madison  6,  Wisconsin 

1 Professor  R.  E.  Langer 
1 Professor  P.  C.  Hammer 

Professor  Wallace  Givens 
Department  of  Mathematics 
Northwestern  University 
Evanston,  Illinois 

Dr.  A.  3.  Householder 

Oak  Ridge  National  Laboratories 

Oak  Ridge,  Tennessee 


1 


9 


t 


DISTRIBUTION  LIST  (Continued) 

Recipient 

Professor  John  Todd 
Department  of  Mathematics 
California  Institute  of  Technology 
Pasadena,  California 

Dr.  Franz  Alt 

Mathematical  Computation  Laboratory 
National  Bureau  of  Standards 
Washington  25,  D.  C. 

Dr.  E.  W.  Cannon 
Applied  Mathematics  Division 
National  Bureau  of  Standards 
Washington  25,  D.  C. 

Miss  Irene  Stegun 

Mathematical  Computation  Laboratory 
National  Bureau  of  Standards 
Washington  25,  D.  C. 

Professor  John  W.  Carr,  III 
Director  Research  Computation  Center 
University  of  North  Carolina 
Chapel  Hill,  North  Carolina 

Professor  J.  B.  Rosser 
Department  of  Mathematics 
Cornell  University 
Ithaca,  New  York 

Professor  Richard  S.  Varga 
Department  of  Mathematics 
Case  Institute  of  Technology 
Cleveland,  Ohio 

Professor  D.  H.  Lehmer 
Department  of  Mathematics 
University  of  California 
Berkeley  4,  California 


